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ABSTRACT 


The  general  conditional  equations  which  govern  the  phase  equilibria 
in  three-component  systems  are  presented.  Using  the  general  conditional 
equations,  a  general  method  has  been  developed  to  precalculate  ths  phase 
equilibria  in  three -component  systems  from  fir st  principle  using  computer 
technique.  The  method  developed  has  been  applied  to  several  model  examples 
and  the  system  Ta-Hf-C.  The  phase  equilibria  in  three-component  systems 
calculated  using  the  simplified  method  as  originally  developed  by  Rudv,  agree 
well  with  those  calculated  by  the  present  method.  The  only  difference  is  in 
the  homogeneous  range  with  respect  to  the  interstitial  component  of  solid  solu¬ 
tions  which  exhibit  large  variation  with  metal  exchange.  This  is  to  be  expected 
in  view  of  the  assumptions  made  in  the  simplified  method. 

In  connection  with  the  phase  diagram  calculation  and  other  problems  of 
the  present  Air  Force  contract,  several  computer  programs  have  been  develop¬ 
ed  which  are  included  in  the  appendix  of  this  report. 
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I. 


INTRODUCTION  AND  SUMMARY 


A.  INTRODUCTION 

Apart  from  theoretical  interest,  the  capability  of  predicting 
phase  equilibria  in  three-component  systems  from  binary  data  is  of  practical 
significance  since  the  cost  of  experimentally  investigating  phase  equilibria 
in  ternary  systems  is  rather  high.  For  the  transition  metalj,  metal2 .carbon 
ternary  systems  (i.e.  the  interstitial  type  of  solid  solutions),  a  method  has 
been  developed  by  Rudy^  to  pred.ct  the  phase  equilibria  in  ternary  systems 
assuming  that  the  intermediate  phases  in  the  meta^-carbon  and  meta^-carbon 
binaries  are  either  of  line -compounds  or  that  both  phases  are  of  equal  depend¬ 
ence  of  the  free  energy  on  the  concentration  coordinate  of  the  interstitial 
component.  In  cases,  where  the  homogeneous  ranges  of  the  solid  solutions 
with  respect  to  the  interstitial  component  changes  drastically  with  metal 
exchange,  the  simplified  method  does  not  predict,  as  to  be  expected  from  the 
assumptions  made,  the  exact  phase  boundaries.  The  purpose  of  the  present 
work  is  to  develop  a  general  method  using  computer  technique  to  predict  the 
phase  equilibria  in  three -component  systems  with  no  assumptions  made  as  in 
the  simplified  method.  This  general  method  will  predict  not  only  the  phase 
equilibria  but  also  the  correct  homogeneous  ranges  of  the  single  phases. 

B.  SUMMARY 

The  conditional  equations  which  govern  the  phase  equilibria  in 
three-component  systems  are  presented.  Using  the  conditional  equations,  a 
method  has  been  developed  using  computer  technique  to  precalculate  the  phase 
equilibria  in  three- component  systems  from  first  principle.  The  method 
developed  has  been  applied  to  model  examples  and  the  system  Ta-Hf-C.  The 
phase  equilibria  in  three -component  systems  calculated  using  the  simplified 
method  originally  developed  by  Rudy,  agree  rather  weil  with  those  calculated 
by  the  present  method.  The  only  difference  ;s  in  the  homogeneous  ranges  of 
solid  solutions  which  exhibit  large  variation  with  metal  exchange.  This  is  to 
be  expected  in  view  of  the  assumptions  made  in  the  simplified  method. 


i 


Several  computer  programs  for  the  phase  equilibrium  calcu¬ 
lations  as  well  as  for  other  problems  of  rhe  present  Air  Force  contract  have 
been  developed.  The  fortran  statements  of  all  the  computer  programs  are 
included  in  the  appendix  at  the  end  of  the  report. 

II .  THERMODYNAMIC  DESCRIPTION  OF  PHASE  EQUILIBRIA  IN 

THREE  COMPONENT  SYSTEMS 

According  to  the  phase  rule,  at  constant  pressure  the  maximum  number 
of  coexisting  phases  in  a  one-component  system  is  two,  two-component  system 
three,  and  three -component  system,  four.  In  addition,  if  one  fixes  tempera¬ 
ture,  the  maximum  number  of  coexisting  phases  in  a  ternary  system  reduces 
to  three.  Consequently,  an  isothermal  section  of  a  ternary  phase  diagram  is 
built  up  of  one-phase,  two-phase,  and  three-phase  equilibria.  On  the  one  hand, 
concentrations  of  the  coexisting  phases  in  a  two-phase  field  determine  the 
phase  boundaries  of  single-phase  fields.  On  the  other  hand,  the  boundaries  of 
a  three-phase  equilibrium  are  the  limiting  tie-lines  of  the  three  adjacent  two- 
phase  equilibria.  We  shall  now  discuss  first  the  thermodynamics  of  two-phase 
equilibria. 


A. 


TWO-PHASE  EQUILIBRIA  IN  THREE-COMPONENT  SYSTEM 


The  Gibbs  free  energy  of  formation  of  a  two-phase  alloy, 

AxByGz  as  shown  in  Figure  1  is  a  linear  combination  of  the  Gibbs  free  energies 
of  formation  of  the  two  coexisting  phases  A^jB^,  G^,  and  AxllB^,,C  Express¬ 
ing  in  mathematical  terms,  we  have 


AG  r  VjAGj  +  v2AG  2  (1) 

where  AG,  AG,,  and  AG?  are  the  Gibbs  free  energies  of  formation  of  the  two- 
phase  alloy  and  ot  the  two  co-existing  single  phases;  and  v,  and  v2  are  the 
relative  amounts  of  the  twe  co-existlng  single  phases.  If  we  take  one  gramatom 
of  alloy  as  our  basis,  we  have  the  following  three  bo-  ndary  conditions: 
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1 


(2) 

(3) 

(4) 


x '  +  v '  +  z 1  -  1 

X"  +  y"  +  Z"  -  1 


Two  additional  boundary  conditions  can  be  obtained  resulting  from  the  con¬ 
servation  of  masses: 


l^X  :  + 

V  X* 

2 

'  X 

v,y'  + 

v2yM 

x  y 

(5) 

(6) 
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Given  AGj  and  AG2  as  a  function  of  composition  at  constant  P  and  T  and  a 
set  of  x(  y,  z  values,  we  would  like  to  calculate  the  values  of  x'.y'.z'  and 
x",y",  and  z".  Since  there  are  eight  unknowns:  vf  ,  v2,x y'»  y",  and  z" 

and  five  equations  (2  to  6),  we  have  only  three  independent  variables.  We  can 
choose  any  three  variables  as  we  wish.  In  the  present  case,  we  shall  choose 
x',y',  and  x".  as  our  ultimate  variables.  By  elimination  of  the  various  vari¬ 
ables  in  equations  (2)  through  (6),  we  obtain  the  following  expressions  for 
Vj,  v2,  z',  y"  and  z"  in  terms  of  x',  y'  and  x": 


and 


X  -  x" 

(?) 

V1  “ 

x'-  x" 

V2 

x  -  x* 
x"-  x* 

(8) 

z'  = 

1  -  x' 

-  y‘ 

(3a) 

v*  •  = 

x'-  x" 

X  -  x"  . 

(9) 

y 

A. 1  -  X 

^  X  -  X1  ^ 

z"  = 

1  -  X" 

-  y” 

(4a) 

In  order  to  solve  for  x^y'  and  x”,  we  need  three  additional  equations  which  must 
be  derived  from  the  condition  that  at  equilibrium  AG  according  to  equation  (1) 
is  a  minimum.  We  can  minimize  AG,  with  the  five  constraints  expressed  by 
equations  (2)  through  (6),  after  the  method  of  Lagrange  as  was  originally  done 
by  Rudy^ .  We  obtain. 
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(15) 


(16) 


(17) 


From  equations  (12),  (13),  (15),  and  (16);  (12),  (14),  (15),  and  (17);  and  (12), 
(11),  (15),  and  (17),  we  obtain  the  following  additional  equations: 


[aAG! 

3ag; 

raAG2 

3AG] 

1  8x' 

3yf 
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AG 


(18) 


(19) 


fa  AG,  3  AG  ]  [BAG,  9  AG] 

r  AG*  = <**-  *">  hr  -  -^-JTfP+  {y'- y,,)  hr-  -  *HTiP  {20) 


From  equations  (18),  (19),  and  (20),  we  can  in  principle  solve  for  x",  x',  and 
y1.  However,  since  the  equations  for  AGj  and  AG2  are  non-linear  with  respect 
to  the  three  independent  variables  x”,  x1,  and  y1,  the  solving  of  the  three 
simultaneous  equations  is  not  very  simple.  Instead,  it  is  more  convenient  to 
numerically  find  the  minimum  value  of  AG  according  to  equation  (1)  by  varying 
the  values  of  x’’,  x',  and  y'  from  0.0  to  1.0.  We  can  take  any  small  incremental 
value  as  we  *ish.  In  general,  for  phase  diagram  calculations,  we  have  found 
a  value  of  0.01  is  sufficient. 

Since  x1  t  y'  s  1,  v.e  need  only  to  consider  the  values  of  x',  y', 
and  x"  in  the  prism  as  shown  in  Figure  2  for  the  computation  of  AG.  Moreover 
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Figure  2.  Schematic  Representation  of  the  Three  Independent 
Variables  x',  y'  and  x"  for  Two-Phase  Equilibria 
in  a  Ternary  System  when  x"  >  x. 


the  values  of  x  always  lie  between  those  of  x1  and  x" .  For  each  value  of  x” 
greater  than  x,  we  compute  AG  for  all  values  of  x1  smaller  than  x,  i.e.  the 
shaded  trapezoidal  region  (in  Figure  Z)  and  for  each  value  of  x1  smaller  than 
x,  we  do  the  similar  computations  for  all  values  of  x'  greater  than  x,  i.e.  the 
shaded  trapezoidal  region  in  Figure  3. 

Since  v  » and  y”  are  all  positive,  rearranging  equation  (6), 

we  have 
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Figure  3.  Schematic  Representation  of  the  Three  Independent 
Variables  x',y‘,  and  x"  for  Two-Phase  Equilibria 
in  a  Ternary  System  when  x"  <  x 


y 


(21) 


Equation  (21)  gives  an  upper  limit  of  y1. 


The  number  of  computations  for  searching  the  minimum  of 
AG  according  to  equation  (1)  can  be  further  reduced  when  we  use  the  additional 
information  provided  by  the  binary  phase  boundaries. 

Before  proceeding  further,  we  must  have  some  mathematical 
expressions  for  AGJ  and  AG2.  As  shown  previously, the  free  energy  of  a 
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single  phase  of  the  interstitial  type  of  solid  solutions  such  as  the  ternary 

-wHidef  &pPh%zmry 

•  — iv  be  adequately  represented  by  the  following  expression, 

AG,  (in  cal/gatonv  alloy)  -  x'AG^  +  y'AGg^  +  |x'ln^_T  +  y  'inj^rj  (22) 

where  x'(y'  and  z'  are  atom  fractions  of  A,  B,  and  C  in  phase  1,  «'  is  the 

interaction  parameter  for  the  solid  solution  (A,B)Cu»  R  is  tKe  universal  gas 

constant,  T  is  the  absolute  temperature,  and  AG.  -  and  AG-,-  are  the  Gibbs 

AGU  BCU 

free  energies  of  formation  of  the  two  respective  binary  phases  in  cal/gatom  A 
or  B.  We  have  a  similar  expression  for  AGZ  which  is  a  function  of  x",  y"  and  z" , 

Often  we  have  the  cases  where  the  solubility  of  the  third  com¬ 
ponent,  let's  say  C,  in  the  solution  (A,B)C^  is  so  small  that  for  practical  pur¬ 
poses,  we  can  take  u  as  zero.  For  such  cases,  z1  =  0,  we  reduce  the  three 
independent  variables  x',  y'  and  x"  to  only  x'  and  x".  Therefore,  the  search 
for  a  minimum  for  AG  according  to  equation  (1)  reduces  from  a  three-dimensional 
to  a  two-dimensional  problem. 

Model  Example  1: 

Let  us  now  take  a  hypothetical  case  where  A  and  B  form  a  series 
of  continuous  solid  solutions  at  high  temperature;  A  and  C  and  B  and  C  form 
two  intermediate  phases  AC  ,  AG  and  BC  ,  BC  ,  respectively  as  shown  in 

r  V  W  V  w  r  7 

Figures  4  and  5.  Moreover,  the  solubility  of  C  in  the  solid  solution  (A,B)  is 
assumed  to  be  small;  and  (A,  B)Cv  and  (A,B)Cw  form  two  series  of  continuous 
3olid  solutions.  The  homogeneous  range  with  respect  to  component  C  is  narrow 
in  the  solid  solution  (A.BJC^  solid  solution.  We  would  like  now  to  compare  the 
tie-line  distributions  in  the  two-phase  regions  (A,  B)-(A,  B)C  and  (A,  B)-(A,  B)G 
calculated  at  2000*K  using  the  method  described  here  with  those  calculated  by 
the  simplified  method^ 

The  Gibbs  free  energies  of  the  four  binary  intermediate  phases; 

AC^,  BG^  ,  AC^f,  and  BC^  may  be  represented  by  the  following  equations  (see 
Figures  4  and  5), 
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ATOMIC  FRACTION  OF  C 


4.  Gibbs  Free  Energies  of  the  Two  Intermediate  Phases 

AC  and  AC  in  the  oinary  A-C. 
v  w  1 


AG 


AC 


AG 


BC 


{23,  000-240,  000  z  +  600,  000  i  } 
J—  (35,  710-418,  882  z  +  1,  102,  162  z*  } 


ATOMIC.  FRACTION  OF 


Figure  5.  Gibbs  Free  Energies  of  the  Two  Intermediate  Phases 

BC  and  BC  in  the  Binary  A-C. 

v  w  7 


AGAC  =  TTI  {121 , 0f> 5-504,  220  z  +  504,  22  0  z*  } 
w 

aCBC  5  {15,  335-83.  555z  +  89,  780  z*  } 

w 

The  Gibbs  free  energies  of  the  ternary  solutions  (A,B)Cv  and  (A,  B)C  m 
take  the  same  form  as  equation  (22).  The  interaction  parameters  for  the 
three  solid  solutions  are  taken  to  be, 


« 


« 


( 


A,  B 

=  5970 

(A,  B)C  v 

=  3980 

(A.B)Cw 

=  3980 

We  shall  first  consider  the  two -phase  *  quilibria  between  the 
solid  solutions  (A,B)  and  (A,  B)Cv»  Since  we  assume  that  the  intermediate 
phases  AC^  and  BC^  are  not  stable  in  the  binary  systems,  the  phase  AC 
then  becomes  stable.  Accordingly,  we  have 


AG  =  5970  x'y'  +  RT  (x1  In  x'  +  y'  In  y') 


AC,  .  *..AGAC  *  y"AGBC  .  * 


|x" 


+  RT  I x"ln  +  y"ln 


_lll 

l-z”J 


(27) 


(28) 


One  can  see  from  Figures  6  and  7  that  the  tie-line  distribution  between  the 
solid  solutions  (A,B)  and  (A  B)Cv  calculated  using  the  method  described  here 
agrees  well  with  that  calculated  using  the  simplified  method.  This  is  to  be 
expected  since  the  range  of  homogeneity  with  respect  to  component  C  in  (A,B)Cv 
is  relatively  narrow.  The  simplified  method  for  calculating  the  two-phase 
equilibria  in  ternary  systems  has  been  discussed  extensively  by  Rudy^.  It 
suffices  to  say  that  the  compositions  of  the  two  co-existing  phases  are  where 
the  concentration  free  energy  gradients  are  the-  same,  i.e. 


dAG, 

“Hx7-- 


dAG2 

dx" 


(29) 


Figure  6.  Phase  Equilibria  Between  Phases  P  and  7  and  7  and  C 
in  the  System  A-B-G  at  2000*K  Calculated  Using  the 
General  Method. 

The  concentration  free  energy  gradient  curves  for  the  solid  solutions  (A,B), 
(A,B)Cv  and  {A,B)Cw  are  shown  in  Figure  8. 

Next,  we  consider  the  two-phase  equilibria  between  ti  t  solid 
solutions  (A,B)  and  (A,B)Cw.  In  this  case,  we  assume  that  the  phases  AC^ 
and  BCv  do  not  appear  in  the  two  respective  binaries.  Since  the  homogeneous 
range  for  the  phase  BC ^  ia  relatively  large  and  the  free  energies  of  AC ^ 
phase  increase  rapidly  at  concentrations  of  C -component  higher  and  lower 


Figure  7.  Phase  Equilibria  Between  Phases  (3  and  -j  and  >  and  C 
in  the  System  A-B-C  at  2000*K  Calculated  Using  the 
Simplified  Method. 


0.5,  the  tie-line  distribution  calculated  (Figure  9)  used  this  simplified  method 
does  not  agree  well  with  that  calculated  (Figure  10)  using  the  general  method 
described  here  ag?in  as  expected.  Moreover,  the  A-B-rich  phase  boundary 
of  the  solid  solution  (AtB)Gw  calculated  using  the  method  described  in  this  paper 
shows  a  curvature  as  expected  while  that  calculated  using  the  simplified 
method  does  net. 

Ii  should  be  pointed  out  here  since  the  solubilities  of  A  and  B 
in  C  are  assumed  to  be  negligible,  the  tie  lines  in  the  two-phase  regions 
(A,B)Gv  and  C  and  (A,B)Cw  and  C  (Figures  6,  7,  9,  and  10)  were  not  calculated 
and  must  be  pointed  toward  the  pure  component. 
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Figure  8.  Concentration  Free  Energy  Gradient  Curves  for 

Solutions  (A ,B),  'A,  B)Cv  and  (A,B)Cw  at  2000*K. 


Figure  9.  Phase  Equilibria  Between  f3  and  6  and  6  and  C  in 
the  System  A-B-C  at  2000*K  Calculated  Using  the 
Simplified  Method. 


B.  THREE-PHASE  EQUILIBRIA  IN  THREE-COMPONENT 
SYSTEM 


In  a  similar  manner  as  for  the  two-phase  equilibria,  the  Gibbs 

free  energy  of  a  three-phase  alloy  A^B^C^  as  shown  in  Figure  11  is  a  linear 

combination  of  the  Gibbs  free  energies  of  formation  of  the  three  co-existing 

phases  A  ,B  ,C  ,,  A  ,,B  ,,C  and  A  ,,.B  , ,,C  MI: 
r  x'  y'  z'  x"  y"  z"  x'"  y"'  z1" 


15 


Figure  10.  Phase  Equilibria  Between  (3  and  6  and  6  and  C 

In  the  System  A-B-C  at  2000*K  Calculated  Using 
the  General  Method. 


AG  -  VjAGj  +  v2AG2  +  v^AGj 

(30) 

with  the  following  constraints: 

Vj  +  w2  +  V,  =1 

(3D 

x'  +  y1  +■  z'  =1 

(32) 

Figure  11.  Three-Phase  Equilibrium  in  a  Ternary  System 


x"  +  y1'  +  z1' 

=  i 

(33) 

x  "  '+  y ' 1 '+  z'" 

=  i 

(34) 

y  X  1  +  y  XM  +  V  X  1  *  1 

1  2  3 

-  X 

(35) 

y^V*  +  y^yn  +  v  y 1  * 1 

=  y 

(36) 

We  have  twelve  unknowns  and  six  equations  (31  through  36).  Therefore,  we 
have  six  independent  variables.  It  is  no  longer  practical  for  us  to  use  the 
similar  method  as  for  the  two-phase  equilibria  to  find  the  combination  of  the 
compositions  of  the  three  co-existing  phases  where  AG  according  to  equation 
(30)  is  a  minimum. 
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Fortunately,  we  can  reduce  the  six  independent  variables  to 
two  by  taking  advantage  of  the  fact  that  the  boundaries  of  a  three-phase 
equilibrium  are  nothing  but  the  limiting  tie-lines  of  the  three  adjacent  two- 
phase  equilibria.  We  shall  first  discuss  the  method  of  seeking  the  composi¬ 
tions  (x‘",  y*'1,  *'*')  of  the  third  co-existing  phase  when  the  compositions 
(x',  y',  z'  and  x",  y",  z")  of  the  other  two  co-existing  phases  are  known.  We  will 
then  discuss  how  to  find  the  compositions  x'.y'.z'  and  x",y",z".  Under  such 
conditions,  instead  of  having  the  six  constraints  as  expressed  by  equations 
(31)through  (36),  we  have  only  four  boundary  conditions:  Equations  (31),  (34), 
(35),  and  (36)  with  six  unknowns,  i.e.  Vj,  v2 ,  v^,  x" ',  y ' 1 '  and  x' 1 ' .  Therefore, 
we  have  ultimately  two  independent  variables.  We  can  choose  any  two  variables 
as  we  wish,  but  from  a  practical  point  of  view,  it  is  convenient  to  choose  the 
two  variables  x"1  and  y'".  Expressing  ^  ,  v2,  v3  and  x'"  in  terms  of  y1”  and 
z  " 1 ;  we  have 

w  (y1-  y)(z'"-  -  (»"-  z)(y"'-  y) 

1  D 

„  -  (y  -  y')(z"'-  z')  -  (z  -  z  ')(y" '  -  y') 

2  '  D 

Vj  ,  i  -  v,  -  v2 
x'"  =  1  -  y ' 1 1  —  Z1" 

where 

d  =  (y"~y ')(*' ' ' -  *’)  -  (z  -  z’Hy'"-  y') 

To  find  the  values  of  y'"  and  z'"  we  again  have  to  use  the  condition  that  at 
equilibrium  AG  according  to  equation  (30)  is  a  minimum.  The  method  for 
seeking  the  values  of  y'"  and  z'"  where  AG  is  a  minimum  is  similar  to  the 
method  used  for  solving  the  two-phase  equilibria  discussed  earlier. 

For  the  three-phase  equilibria  resulting  from  a  miscibility  gap 
in  one  of  the  solid  solutions,  the  compositions  of  the  miscibility  gap  of  this 
solid  solution  are  well-defined.  In  this  case,  we  can  solve  for  the  composi¬ 
tion  of  third  co-existing  phase  using  the  method  described  here. 


(37) 

(38) 

(31a) 

(34a) 

(39) 
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Model  Example  2. 

We  shall  use  the  same  data  as  In  Model  Example  1  to  illustrate 
the  application  of  the  method  described  here.  When  the  temperature  is  lower 
from  2000*K  to  1250®K,  a  miscibility  gap  results  from  the  solid  solution  (A,B). 
The  compositions  of  the  miscibility  gap  at  1250’K  are  x'  -  0.83,  y'  -  0.17 
and  x"  =  0.  17  and  y"  =  0.83.  Again,  as  shown  in  Figures  12  and  13,  the  phase 
equilibria  calculated  using  the  simplified  method  agree  well  with  those  cal¬ 
culated  using  the  general  method  with  the  exception  that  homogeneous  range  of 
the  solid  solution  (A,B)Cv  with  respect  to  component  C  is  r>a  •rower  The  free 
energy-concentration-g  adient  curves  used  to  determine  the  phase  equilibria 
at  12bO®K  as  shown  in  Figure  12,  are  displaced  in  Figure  14. 

Referring  to  Figure  lb,  the  method  for  seeking  the  composition 
of  the  three  co-existing  phases  P-7-&  will  be  discussed  using  the  following 
conditions: 


at  Xj.yj.z,: 

AG  (2$ 

between  P-6) 

A  A 

AG(2q> 

AG(3(j> 

between  3  - 
between  p-7-6) 

at  xz*  V  V 

AG  (3p 

between  p-6-7) 

< 

< 

AG(2<$> 
AG  (2$ 

between  p-6) 
between  p  -7) 

at  xj*  y3.S 

AG  (2$ 

between  p-7) 

< 

< 

AG(2y  between  P-6) 
AG(3V  between  p-7-6) 

In  reality,  we  don't  have  any  idea  about  the  compositions  of  the  three  co-existing 
phases  except  we  know  there  is  a  three-phase  equilibrium  p-7-6  since  the 
corresponding  phase  in  the  A-C  binary  is  not  stable.  Thus,  we  can  start,  let's 
say,  with  a  composition  Xj  ,  yJt  Zj  and  ask  ourselves  what  possible  three-phase 
field  it  is  in.  We  know  one  of  the  phase  boundaries  of  this  three-phase  equi¬ 
librium  |3- 7-6  is  the  tie-line  between  (3-6  and  furthermore,  this  tie-line  must 
lie  at  the  left  of  the  composition  x,  *  Yj  >  •  Thus,  we  can  start  with  various 

compositions  at  the  left  of  the  composition  .  y, .  z,  and  first  find  the  correspond¬ 
ing  compositions  of  the  two  co-existing  phases  p  and  6.  Knowing  the  various 
values  of  x'.y'.z'  (p -phase)  and  x'  ,  y",  z"  (6  -phase),  one  can  calculate  the 
corresponding  values  of  x'",  y"',  z  "  (7-phase)  and  the  corresponding  values 
of  AG  (3$  between  p-6-7)  .  One  then  compares  the  lowest  value  among  all  the 
AG  values  calculated  with  those  of  AG  (2$  between  p-5)  and  AG(2$  between  p-->) 
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Figure  12.  The  Phase  Equilibria  Calculated  Using  the  Simplified 
Method  in  the  System  A-3-C  at  1251 *K  When  the  Two 
Binary  Phases  ACw  and  BCw  are  Unstable. 

using  the  method  described  earlier.  The  stable  equilibrium  is  defined  by 
the  lowest  value  of  the  Gibbs  free  energy.  If  AG  (2$  between  P-6)  is  the 
lowest,  thu  three-phase  equilibrium  must  lie  at  the  right  of  the  chosen  values 
of  X! '  y i  <  zj  ■  H  on  the  other  hand,  the  value  of  AG  (2<$>  between  P-y)  is  the  lowest, 
then  the  three-phase  equilibrium  must  lie  to  the  left  of  the  chosen  values  of 
Vy  .z  .  If  the  value  of  AG  (3$  between  p-y-6)  is  the  lowest,  than  we  have 
found  the  three-phase  equilibrium.  A  simple  and  rapid  iteration  process  can 
be  prepared  to  solve  the  compositions  of  the  three  co-existing  phases. 
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Figure  13.  The  Phase  Equilibria  Calculated  Using  the  General 
Method  in  the  System  A-B-C  at  1250*K  When  the  Two 
Intermediate  Phases  AC  and  BC  are  Unstable. 

w  w 


Model  Example  3: 

Using  the  same  data  as  in  Model  Example  1  and  assuming  all 
the  three  phase  BCv>  AC w  and  BCw  appear  in  the  corresponding  two 
binaries  A-C  and  B-C*  we  shall  use  the  method  described  here  to  find  the 
three-phase  field  £->-6  and  the  tie-lines  in  the  three  two-phase  fields:  j3-6, 

P-7  and  7-6.  The  phase  diagrams  calculated  at  2000*K  using  the  general 
method  as  shown  in  Figure  16  agrees  well  with  that  calculated  using  the  simpli¬ 
fied  method  (Figure  i  7} .  The  obvicus  difference  between  the  two  calculated 
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2.0 


1.  (A-B)  Solid  Solution 

2.  (A-B)  Cv  Solid  Solution 


Ficure  14.  The  Concentration-Free -Energy-Gradient  Curves 

for  Solutions  (A>B)  and  (A,B)Cv  at  1250  K. 


A 


Figure  17.  Phase  Equilibria  Calculated  Using  the  Simplified 
Method  in  the  System  A-B-C  at  2000*K  When  the 
Three  Phases  AC,.  BC  and  BC  are  Stable  in 
the  Binaries.  v  w 


phase  diagrams  is  the  (A-B)  rich  homogeneous  range  of  the  solid  solution 
(A-B)GW  with  respect  to  C -component.  The  homogeneous  range  calculated 
using  the  general  method  shows  the  curvature  while  that  calculated  using  the 
simplified  method  does  not. 

Although  the  technique  of  fixing  the  compositions  of  the  co¬ 
existing  phases  for  a  three-phase  equilibrium  using  the  simplified  method 
has  been  discussed  extensively  previously,  it  is  worthwhile  to  point  out  that 
in  addition  to  the  condition  at  equilibrium,  the  free-energy-concentration 
gradients,  i.e. 
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(A,  B)C . 


ax 
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must  ill  be  necessarily  equal,  we  need  the  additional  stability  condition  in 
order  to  fix  the  compositions  of  the  three  phases.  Expressing  in  thermo¬ 
dynamic  terms,  the  stability  condition  simply  says  that  at  equilibrium  the 
integral  Gibbs  free  energy  of  decomposing  an  alloy  (A,B)Cv  into  two  other  alloys 
(A,B)Cu  and  (A,B)CW  must  be  zero.  For  the  modal  example  described  here, 
the  total  Gibbs  free  energy  of  decomposition  of  the  alloy  (A,  B)C0  ri5  to  the 
alloys  (A,B)C096]  and  (A,B)  according  to  the  following  reaction  is  shown  in 
Figure  18. 


(A,B)C0iIj5  <ss.> -* 


0.244  (A.BJCfcttj  <ss>  +  0.755(A,B)<S£> 


III.  APPLICATION  OF  THE  METHOD  DKVELQFED  TO  THE  SYSTEM 
Ta-Hf-C 


As  already  discussed  earlier,  in  order  to  precalculate  the  phase  equi¬ 
libria  in  a  ternary  system,  it  is  necessary  to  know  the  free  energies  of  the 
binary  phases  as  a  function  of  composition  and  the  solution  behaviors  of  the 

ternary  phases.  In  the  binary  system  Hf-C,  one  intermediate  phase,  hafnium 

(4) 

monocarbide  HfC  (Bl-type),  forms  with  a  wide  range  of  homogeneity1  '.  The 
a-Hf  terminal  phase  is  stabilized  by  the  addition  of  carbon  atoms  to  lusher  tem¬ 
peratures  However,  since  little  is  known  about  the  free  energy  of  this  phase 
as  a  function  of  composition,  we  shall  at  the  present  ignore  the  stabilization  of 
the  a-phase  to  high  temperature s  and  assume  that  metal  and  carbon  components 
at  the  metal-rich  boundaries  of  the  monocarbide  phase  are  in  equilibrium  with 
those  in  the  |3-Pf  terminal  phase  at  high  temperatures. 

In  addition  to  the  rnonocarbide  phase  (Bl-type),  a  second  intermediate 
chase.  Ta2G  forms  in  the  system  tantalum-cerbor..  Tantalum  subcarbide,  TazC, 
is  polymorphic.  The  low -temper  a  tore  form,  a-Ta  G,  has  a  hexagonal  close- 
packed  arrangement  of  the  metal  atoms.  The  distribution  of  the  carbon  atoms 
among  the  interstitial  sites  is  not  known  but  undoubtedly  dependent  on  temperature. 
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The  high-temperature  form  p-Ta^C,  probably  has  tne  same  arrangement  of 
metal  atoms®^  but  differs  in  the  degree  of  disorder  in  the  carbon  sublattice 
from  the  a -modification.  Again,  for  the  lack  of  pertinent  thermodynamic 
data  with  respect  to  the  two  modifications  of  Ta2C,  we  shall  assume  a-Ta^C 
and  jl-Ta^C  to  be  of  one  phase  for  the  calculation  at  the  present. 

Tne  Gibbs  free  energies  of  the  monocarbide  and  subcarbide  phases 
will  be  represented  by  a  sum  of  two  contributions:  the  thermal  and  the  con¬ 
figurational  contributions  to  the  total  free  energies  of  these  phases,  i.e. 

AG  =  AGth  +  AGC  ^  (39) 


Since  the  free  energies  of  formation^  for  both  of  the  two  monocarbide  phases 
at  the  stoichiometric  composition  have  been  determined,  we  shall  represent 
the  free  energies  of  formation  due  to  thermal  vibration  for  compositions  less 
than  x^.  -  0.  5  by, 

AG**1  (B1  phase)  -  Z(  l-z)  AG^  +  (l-2z)(a  +  bz  +  Cz2)  (40) 

where  AGf  is  the  Gibbs  free  energy  of  formation  of  the  stoichiometric  mono¬ 
carbide  phase  in  cal  gatom  alloy,  z  is  the  atom  fraction  of  the  carbon  com¬ 
ponent,  and  a,  b  and  c  are  parameters  needed  to  be  determined.  When  d  and  c 
are  set  equal  to  zero,  we  have  the  Wagner -Schottky^’  8’  ^  vacancy  model  with 
'a'  being  the  Gibbs  free  energy  of  forming  a  carbon  vacancy  on  the  carbon  sub¬ 
lattice.  For  the  free  energy  contribution  due  to  configurational  mixing  of  carbon 
atoms  among  the  various  lattice  sites,  we  assume  ideal  entropy  of  mixing  which 
is  reasonable  at  high  temperatures.  Accordingly,  we  have 


AGconf  -  RT 


[zln  ihi 


l 


+  (l-z)  In 


1  -2z* 
l-z 


(41) 


For  the  subcarbide  phase,  the  thermal  contribution  to  the  total  free 
energy  will  be  represented  by  a  3rd  degree  polynominal,  i.e. 


AG  -a'  +  b'z+c'z  +d'z5 


(40b) 
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For  the  configurational  free  energy,  we  shall  use  the  model  developed  by 
Rudy^*^.  In  this  model,  one  assumes  that  the  two  interstitial  sites  occupied 
by  the  carbon  atoms  are  energetically  different  by  an  amount  of  AE.  Express¬ 
ing  in  terms  of  two  gram  atoms  of  metal  (Me  C  1.  the  configurational  free 
energy  as  derived  earlier'  '  is 


AGMe2Cw  *  ZBAE  +  Kln  V  ll-ZA>  ln 


+  z0ln  z0  +  (1-Zg)  In  0-zB)]  RT  (42) 


where  AE  is  the  energy  difference  between  the  two  interstitial  sites  in  the 
cal/2 -gram  atom  metal,  is  the  mole  fraction  of  carbon  atoms  on  the  A-sites, 
(1-*A>  is  the  mole  fraction  of  vacant  A-sites,  is  the  mole  fraction  of  the 
carbon  atoms  on  the  energetically  unfavorable  B-sites  and  (1-z  )  is  the  mole 
fraction  of  the  vacant  B-sites.  It  is  understood  that 


ZA  +  ZB 


(43) 


However,  since  it  is  more  convenient  to  discuss  the  free  energies  of  solid  solu¬ 
tions  in  terms  of  one  gram  atom  alloy,  we  shall  divide  equation  (42)  by  (2  +  w) 
and  change  the  variable  from  w  to  z  (the  atom  fraction  of  carbon  atoms  in  the 
solution).  Accordingly,  we  have 


4GMe,f'  C  ZBAE  +[ZA  111  ZA  +  (1-ZA}  ln  (1"ZA>  + 

1  -z  z  v 

4  ZB  ln  ZB  *  (1-ZB)  ^  <1"ZB)]  RT]  (44) 

Using  the  well-known  thermodynamic  relationships  relating  the  partial  molar 
and  integral  quantities,  we  have. 
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r-  +  ~  [ZAln  ZA  +  {1"ZA>  ln  (1'ZA>  + 
+  «B1  n  *B  +  <l-zB)  ln  (l-zB)J 


, -  AE  s  +  RT  r  ln  +  RT  s  ln 

i-z  I 


h.] 

^bJ 


conf 

AGc 


AG  +  (1-z) 


dAG 


where 


dz 


z  z 

^  AE  s  +  r  RT  ln  +  s  RT  ln  ~ 

1“Z ,  1“ zr 


(l-zA)Z(l  +  u)2 


rr:  ,2-aE7rTt 

(i-zA)  (1+")  +e 


(l-Zg)2  (1  +  v)2 


AE/RT 


(1-zA)  (1TU)  +e 


AE 
'  RT 


AE 


RT 


1~ 


B 


AE 

AE 

l-3z  ,  „  RT 

IJ+Z1 

n-3z  RT 

1-z  e 

1 1  -  zl 

l  1-Z  T  e 

1+zl "  I  8z  l 

tttJ  1fir)e  1 


AE  RT, 


H  t*l 


(52) 


If  we  set  AE  -  0,  z*  -  z_  since  there  is  no  difference  between  the  carbon 
A  B 

atoms  on  the  A-sites  and  B -sites.  From  equation  (52)  one  obtains. 


z 


A 


z 

1-z 


(53) 


Substituting  for  z  and  z_  according  to  equation  (53),  equation  (44)  reduced 
to  equation  (41). 

Before  we  can  proceed  any  further,  we  must  know  the  values  of  a,b,  c, 
(equation  40a)  and  of  a'jb'jC1  (equation  40b)  for  the  monocarbide  and  subcarbide 
phases  in  addition  to  the  values  of  AE  for  the  two  subcarbide  phases.  These 
parameters  were  determined  by  using  the  following  conditions. 

a.  Values  of  AG,  for  TaC  and  HfG  are  those  selected  by 

Chang(6). 


b.  The  free. energy  decomposition  of  TaCg^  to  0,633 
TaCo.79  +  0.367  Ta  is  taken  to  be  2300  cal  at  2273®K^. 

c.  The  free  energy  difference  between  HfC0>5  and  TaCe>5 
is  taken  to  be  2250  cal  at  2273*K(o). 

d.  The  phase  boundaries  of  the  subcarbide  and  monocarbide 

(4) 

phases  are  those  determined  by  Rudy'  and  Rudy  and  Harmon'  ' . 

e.  The  partial  molar  quantities  in  the  two-phase  regions: 
metal-subcarbide,  subcarbide -monocarbide  and  monocarbide -graphite  must 

be  all  the  same. 

With  these  five  conditions,  we  have  obtained  the  following  results  for 
our  computation: 
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HfG -phase  (Bi-type) 


2AG,  -  -52,  350  +  2.  08  T  298. 15*-2073*K 

G  0.5 

-  -55,630  +  3.66  T  2073*-249i*K 

a  -  2  50,  1 54;  b-  -968,  150;  c  =  1,  113,  170 

Hf2G  -phase 

a'  -  72,445;  b'  -  -466,  000;  c'  -  650,  000;  d'  =  0;  AE^O 
TaC -phase  (Bl-type) 

AG,  =  -35.  335-1.  7949  log  T  +  6.4757  T 

r»  0.5 

a  -  363,  950;  b  =  1,  527,  083;  1,  770,  833 

Ta?C  -phase 

a1  =  129,848;  b’  =  1 ,  336,  1 1  7;  c’  ^  -4,  7:  6,  667;  d'  =  5,  373,  333 

AE  -  4000  cal. 

The  phase  equilibria  in  the  system  Ta-Hf-G  at  2273*K  and  1273*K  are 
calculated  and  presented  in  Figures  19  and  20.  The  phase  diagrams  calculated 
using  the  simplified  method  are  shown  in  Figures  21  and  22  for  comparison.  As 
might  be  expected,  the  marked  difference  between  the  calculated  phase  diagrams 
by  the  two  different  methods  is  the  homogeneous  range  of  the  monocarbide  phase. 
The  concentration-free -energy-gradient  curves  for  the  metal,  subcarbide  and 
monocarbide  phases  as  well  as  the  integral  free  energy  of  decomposition  of  the 
'•subcarbide  phase  into  the  metal  and  monocarbide  phases  at  2273*K  used  to 
establish  the  phase  equilibria  are  presented  in  Figures  23  and  24.  Similar 
curves  used  to  determine  the  phase  equilibria  at  1773*K  are  not  included  in 
this  report. 
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Figure  19.  Calculated  Phase  Equilibria  in  the  System  Ta-Hf-C 
Using  the  General  Method  at  2273  *K. 
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Figure  22.  Calculated  Phase  Equilibria  in  the  System  Ta-Hf-C 
Using  the  Simplified  Method  at  1773*K. 


L 


IN  CAL 


Figure  2-5.  Integral  Free  Energy  of  Decomposition  of  the  Subcarbide 
phase  into  the  Metal  and  Monocarbide  Phases  at  2273’K. 
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IV. 


DISCUSSION 


The  method  described  in  this  rqsort  has  demonstrated  that  the  phase 

equilibria  in  three -component  systems  ca  .  be  adequately  predicted  from 

binary  data  based  on  the  first  principle  provided  that  the  Gibbs  free  energies 

of  formation  of  all  the  binary  phases  are  accurately  known  as  a  function  of 

comDosition  and  temperature.  Moreover,  one  must  also  know  the  solution 

(2.  3) 

behaviors  of  the  ternary  allots.  It  has  be  ?n  demonstrated  previously'  ’  ,  that 

for  the  interstitial  type  of  solid  solution  sue  h  as  the  transition  metal  carbides, 
the  regular  solution  theory  with  proper  values  of  the  interaction  parameters, 
is  adequate  in  most  cases. 

The  general  method  also  demonstrates  that  the  simplified  method 
originally  developed  by  Rudy  for  predicting  phase  equilibria  in  three  -  component 
systems  is  sufficiently  accurate.  In  cases  where  the  homogeneous  ranges  of 
the  ternary  solid  solutions  with  respect  to  the  interstitial  component  such  as 
carbon  are  small,  the  simplified  method  would  predict  the  identical  results 
as  the  general  method.  Only,  when  the  homogeneous  ranges  of  a  ternary  solid 
solution  with  respec  .  to  the  interstitial  element  change  drastic  ally  with  the  metal 
exchange,  wil)  there  be  a  difference  between  the  phase  diagrams  calculated 
using  the  two  different  methods.  In  practice,  the  phase  equilibria  in  a  ternary 
system  predicted  by  the  simplified  method  are  sufficiei.  :iy  accurate  in  most 
cases . 
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APPENDIX  I 


COMPUTER  PROGRAMS  FOR  TERNARY  PHASE 
DIAGRAM  CALCULATIONS  (Ray  Marler) 


! 


S IBFTC  <-l  FULIST 

DIMENSION  X(5P).Y(50).A3(3).A4(3).TITLE(13) 

NAMFL IST/I NPUT/NOXYS* I  PR  I NT • CA SE .DUMP , XL I M ♦ R » T »FPSP ,  EPSPP . A3* A4. 

1  X  .  Y  .  K  1  .  K  2 

NAMEL I ST/DTADMP/L I NE  * X • J *XP » YP . XPP . YPP , Z PP ♦ DG * DG3 .DG4.GNU1 »GNU2  * 
1FXP.FYP.DG1.C1.C2.FC1.FC2.DG2 

1  FORMAT ( 1H139X53HA  ERO  JET-GENERAL  CORPORA!  I  1 

10  N/1H039X53HC  OMPuTlNo  SCIENCES  DIVISION 
2/1H05OX31HS  ACRAMENTO  PLAN  T//1H057X1 7HP  R  0  G  R  A  M 

3  K / 1 H02  7X  7  7HT  WO-PHASE  EQUILIBRIUM  IN  TE 
4RNARY  SYSTE  M/16H0C  U  S  T  0  M  E  R  97X19HP  R  O  G  R  A  M  M 
5  E  R/12H0Y.  A.  CHANG  109X1 1HR.H.  MARLER//) 

READf 5.17> T ITLF  ; 

2  READ  (5. INPUT) 

WRITF(6.1 ) 

WRITE(6,18) TITLF 
WRITE(6. INPUT  i 
LlM  =  XL  I M  +  .005 
NDUYtPl=  DUMP 

NDUMP2  =  INTI DUMP*XLlM+. 005 )  -  NDUMP1*LIM 

RT  =  R*T 

WRITEI6.1) 

WRITFI6.18) TITLF 
DO  14  !  =  1.  NOXYS 
WR  I  TE ( 6  » 16 ) X  I  I)  ,Y(I! 

PASS  =  1. 

L  =  1 
LINF  =  1 

L  I  MX  =  XL  I M  *  XII)  +  .005 
BLIMX  =  LIMX 
JM  AX  =  LlM  -  L IMX 

DO  12  K  =  K 1  *  K  2  j 

CAY  =  K 

XPP  =  CAY/XL  I M 
DO  12  J  =  1.  UMAX 
XJ  =  J 

XP  =  (BLIMX  +  XJ1/XLIM 

ZPP  =  (XP  -  XPP ) *< 1 .-X ( I  I -Y (  I  )  ) / ( XP  -  X(I)) 

IF(ZPP.GE.1..0R.ZPP.LE.C..0R.XP.FQ.XPP.0R.XP.EQ.X< I ) ) GO  TO  12 

DG3  =  ( ( A3 ( 3) *ZPP  +  A3 ( 2 )  l*ZPP+A3( 1) ) /( 1 .-ZPP )  j 

DG4  =  ( ( A4 ( 3 ) *ZPP  +  A4 ( 2 1  )*ZPP+A4(  1 )  ) / ! 1 , -ZPP  ) 

GNU1  =  ( X(  I )-XPP ) / <  XP-XPP ) 

GNU2  =  <  XP-X ( I ) ) / ( XP-XPP) 

I F ( GNUl *  LT  .0. .OR. GNU  1 .GT. 1 . .OR .GNU 2 -LT . 0« . OR.GNU2  «GT  *  1 . ) GO  TO  12 

YP  =  l.-XP 

YPP  =  1.  -XPP  -  ZPP 

call  XLNX( XP.FXP.S12 ) 

CALL  XLNX! YP.FYP.S12 )  j 

DG1  =  <FXP+FYP)*RT+XP*YP*EPSP 
Cl  =  XPP/I l.-ZPP) 

C  2  =  YPP/ [ l.-ZPP) 

CALL  XLNX( C1.FC1 .512) 

CALL  XLNX( C2 . FC? »S12 ) 

DG2= ( ( FC1  +  FC2  )  *RT-*-C  1  *C2*EPSPP  +  C  1  *DG3  +  C  2»DG4 )*(  l.-ZPP ) 

DG  =GNU1*DG1+GNU2*DG2 
IFIDG.GE.O. )GO  TO  12 
IF (PASS.NF . 1. )GO  T0  5 
PASS  =  0. 

GO  TO  6 
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5  IF  (06  ~  OGOlO)6*8,8 

6  DGOLH  «  06 
XPS  »  XP 
YPS  *  YP 
XPPS-  XPP 
YPP$«  YPP 

2PPS»  ZPP 
LSAVE  *  LINE 

8  IF(NDUMP1.NE.K.0R.NDUMP2.NE.J)G0  TO  9 
WRITE  (6.DTADMP) 

GO  TO  11 

9  IF< IPRINT.FO.O)WRITEI6.10)LINE , XP , YP ,XPP . YPP ,ZPP * DG 

10  FORMATdH  15X14, 6E16. 8) 

11  LINF  -  LINE  +  1 

12  CONTINUE 

WRITE (6 ,10 )LSAVF»XPS. YPS .XPPS .YPPS * ZPP S » DGOLD 

14  CONTINUE 

CASE  *  CASE  -  1. 

IF i CASE.GT  *0* ) GO  TO  2 

15  STOP 

16  FORMAT ( 1H0 19X  4KX  *  E 1 5 . 8 . 10X , 4HY  =  E15.8/1H  1 5X4HL 1 NE7X 2HXP 1 4X 
12HYP1 3X3HXPP1 3X3HYPP 13X3HZPP12X5HDEL  G) 

17  FORMAT ( 1 3A6 ) 

18  FORMAT ( 1H027X13A6 ) 

END 
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I 


% !RCTC  <-i  -a  Fut ! CT 

n  1  t  ON  Xl50).YI50!.TiTLE(131.VAI5).VR(5).RA(5>.R'?<5> 

NAMELIST  /I  NP'JT  /.NOXYS  »  I  PR  I  NT,  CASE  ,DUMP,XLlM,R,T  ,EPSP,EPSPP.VA.VB. 

1  X  ,  Y,<1  ,<?  >D3  iDA  .PA  ,P9  .DEACV.DEBCV  ,DCAC'J«CG9C'J 

NAMEL I ST/DTADMP/K , j .XP.YP ,XPP,yPP,ZPP,OG»DGACV *DGBCV»  GNUl *GNU2  » 
1FXP.FYP  ,DG1  .Cl  ,C2  *FC1  .FC2.D62 

1  FORMAT ( 1H139X53HA  EROjET  -  GENERAL  CORPORATI 
10  N/1H039X53HC  OMPUTING  SCIENCES  DIVISION 
2/1H050X31HS  ACRAMENTO  PLAN  T / / 1 H05 7X 1 7 HP  R  0  G  R  A  ,Y 

3  K  / 1 H027X  77HT  w  0  -  D  w  A  S  E  EQUILIBRIUM  IN  T  E 

4  R  N  A  R  Y  S  Y  S  T  r  M/16H0C  U  S  T  0  M  E  R  97X19HP  R  O  G  R  A  M  ** 

5  E  R/12H0Y.  A.  CHANGICSX 1 i HR .H.  MARlER//) 

REA0<  5,171 TITLF 

2  R  F  AO  (5,  INPUT) 

WRITEI6.1 ) 

WRITE (6,16 ITITLF 
WR ITE i 6 . I NPUT  ) 

LlM  =  XlIM  +  .OP 5 
ND'JMDJ=  DUvP 

NDUMP2  =  INTIDUMp*XLlM+.005  !  -  ND'JMP1*LIM 
flag  =  nbumpi 
RT  =  R*T 
WR  I  TE ( 6  » 1 ) 

WRITE(6,18)TITLF 
10  14  I  =  1,  NOXYS 
WRITE (6  *16 ) X (  I  5  ,Yi I  ) 

PASF  =  1. 

L  I  MX  =  XLIM  s  X ( n  +  .005 
BL I mX  =  LlMX 

JM  AX  =  LlM  -  L'l  mx 
DO  12  K  =  X  1  *  '<  2 
CAY  =  K 

XPP  =  C  A  Y  /  X  L  I ,v' 

DO  12  J  =  1,  JM  A  X 

XJ  =  J 

XP  =  ( R  l  I M  X  +  XJ1/XLIM 

IFIXP.EQ.XPP.OR.XP.FQ.XI ! 1  ICO  to  1? 

zpp  =  (  xp  -  xpp ) * ( i ,-x ( I )-v ( i > i /  ;  xp  -  x (  !  i ; 

: F ( ZPP.FO. .5 ' ZPP  =  .499 

IFIZDP. CT.D4.0P.ZPP.LT. 03)00  TO  12 

call  DLGC ( DGAC7.R . T , Z =P , V  A , R A , DE AC V , FL AG  * S 1 2 ! 

CALL  DLGC  (  DGP.CF  ,R  ,  T  *  ZPP  ,  VP  ,RB  ,  DERC  V  *CL  AG  ,  i  1 2  ! 

GNUl  =  ( X ( I )-XPP )/ < XP-XPP) 

GNU2  =  ( XP-X ( I) t /( XP-XPP ' 

I F ( GNU1.LT.0. .OR. GNUl -GT , 1 . .OR .  G,Nu2 . L  T • C • .OR ,G\U2*GT. i . )G0  TC  12 

yd  =  l.-xP 

you  =  1.  —xpp  -  ZF( 

'ALL  XL‘.'XCXP,FXP.*1?’ 

CALL  XL‘!X(YP,FYP,«12) 

DC-1  -  (  FXP  +  FYD  1 -»P  T  +  XP*YP*E  P  =  P  +  0GAC  U*  XP-'0G°ri'*  YP 

r  i  =  X  p  d  /  (  1  .  -  ?  P  p  1 
0  7  =  YPP/(  1  .-ZPP  ) 

OAll  XLNX ( c 1  .FC1 .* 1 2 1 

r  L L  XLNX  !  Cl  tf" .  r  ,  t  li  i 

rr-2=I(-/'l*cC7)*PT  ♦  *E;:pP  +  2  *D0iP  C  V  +  C  1  *00  AC  V  1  5  (  1  »  -2  P  " •'  I 

"■G  =Gnl:  ;  ■  *'*■  ■  7  *  r  : 

i  f  ( cr.  '• .  >  ■■■■  - 

;  F  (  n  c  <;  .  r  .  -L  .  \  r ■  ~  ■- 
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GO  TO  6 

5  IF  106  -  POOL'-' 16.8*8 
A  OGOLO  »  OG 
XP$  ■  XP 

yps  ■  yp 

XPPS-  XPP 
YPP$«  ypp 
ZPPS*  ZPP 

8  IF(nDUHP1.NE.<.OR.NOUVP2.NE.JKjO  to  y 
WRITE  (6,OT*Pmpi 

FLAG  *  0 
GO  TO  1? 

9  I  F  {  I  PRINT,  FO.O  >WPf  Tt  (  6  .  i-1 )  x P  ,  v  p  ,  xPP  .  YPP  ,  2PP  ,  p.G 

10  FORmatUh  19XSF1 A.8.F16.4) 

12  CONT I NUF 

WR  I  TF  t  6  *  1  0  !  XPS»YPS,XPPS,ypp'  ,  ^opc,n(30L0 

14  continue 

case  -  CASE  -  1. 

I FICASF.GT.O. >00  to  2 

15  STOP 

16  FORMAT < IH023X4HX  =  F10.8.5X.4HY  =  FlO.fl/lH  26X2mxp;ax 
12HYP13X3HXPP 1 3X3HYPP 13X3HZPP12X5HDEL  G) 

17  FORMAT ( 13A61 

18  FORMAT ( 1H027X13A6 1 
ENO 
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tinrr^  «-  _•>  fijl  I  ST 

DJvfNSION  X  <  50  ) »  Y ( 50  I »Z ( 50 ) »  T l Tl  E ( 13) ,CDG3 ( 3 > ,CDG4<  3  I *CDG5  t  3 )  . 

■t  roG6  (  1  i 

NAMFl IST/INPUT/  a.B.C.D.R,T.DELTA,EPSP,EPSPP»CDG3 .CDG4, 

Kr»G*> .  r  r>r,6  •  x ,  y  .  ? 

1  FORMAT ; 1H1 39X53HA  EROJET-GE  NEPAL  CORPORATI 
10  N/1H039X53HC  OvPurlNG  SCIENCES  DIVISION 
2/1H050X31HS  ACRAwENTO  PLAN  T / / 1 H057X 1 7HP  R  0  6  R  A  M 
•>  </lu019x  9  3HC.  CNC°AL  TwO-DHASE  E  Q  U  I  L  I  3  R  I 

4  U  «  in  TERNARY  SYSTE  M/16H0C  U  S  T  C  M  E  R  97X 

5 1 9hP  PC  GRAmm£  R/12H0Y  A.  CH ANG 1 09X 1 1 MR ; H.  MARLER//) 

R  F  AD  {  5  »  ’  1TITLE  <  10 
J  FORMAT (  T>  A*  ,  I  ?  1 
7  RPA^  (3, INPUT) 

WRITF  (6.1) 

WR I TE ( 6 .4 ) TITLE. ID 
4  FORMAT ( 1HQ26X1 3Afe  ,25X  I  2  ) 

WRITE  (6. INPUT) 

WRITE  16.1  ) 

WRITE(6.16)TITlE 

16  FORMAT ( 1H026X13A6//  7X1HX1 2X1HY1 2X1  HZ  1 1 X2HXP1 1 X2HYP 1 1X2HZP11 X3HXPP 
111X3HYPP1 1X3HZPP  9X7HDELTA  G/ ) 

QT  =  P  *  T 

no  14  i  =  1,60 

I F ( X ( I ) .FG.O. . AND.YI  1  )  .FQ.O.. AND  .Z (  I) .EO.O. )G0  TO  16 
tDA5S  =  1 

17  IF< A.EQ.X(  I  ) ) A  r  X (  I ) *DELT  A/ 10 . 

XPP  -  A 

6  IF ( XPP  -  X (  I  )  )  6.13.7 

6  XF  =  X (  I  )  +  DELTA 

X°MAX  =  1. 

GO  TO  8 

7  XP  =  r'. 

X  Pw A  X  =  X (  I  ) 

8  ZP  =  r 

12  vp  =  l .  - x  d - Z  p 

I  r  I  XP,  E"'.  XPP  !  r.r>  jn  11 
GNU!  = ( X ( I ) -XPP ) / ( XP-XOP > 

I  F  [  GN' :  1  » GF  •  1  *  •  C7  , C-NU 1  . L  T  . 0  •  )  GO  TO  11 
GNU?  =  1 •  -  GNU 1 
Y  PP  = ( Y ( I) -GNUl *YP ) /6NU2 
Z°P  *1.-XPP-YPP 

IF(ZPP.GE.1..0R.ZPP.LT.0.)G0  TO  11 
call  CUBIC(CDG3.ZPP.DG3»$11) 
call  Ctd  I  C  <  C0C-4,ZP°  .DC-4, 51 1  ) 

Li  1  =  xpp/  (  1 .  -ZPP  ) 

7  =  YPP/(  1  ,-Zpp  ) 

CALL  X L N X ( U I .FUl  .ill  1 
CALL  XL.NX(U?.FU?,TU  ) 

DG?=  (  (  Ful  +FU7  )  *RT+U1  *’ :i?*FPSPP  +U?*DG4  +  lil  *irO  )  *  (  1  . -Z  PP  ) 

CALL  CUpIC (CPG5.ZP.DG5.1 11 ) 

CALL  CUEICICOGt  ,  Z  0  »  ^  C-  6  ,  S  1 1  ! 

U 1 =  XP/(1.-ZP) 

{.'2-  YP/d.-ZD) 

call  XLNX ( Ul ,  n.I  .Til) 

CALL  XLNX  (  U2,  F”.'  ,  i  1  1  ) 

70  l  -  (  (  Fiji  *c  JZ  )  *RT  +  Ul*U2*EP6P  *J?»Dr'6  *  ; ’  l  * ~  c  6  !  *  (  1  .  - Z P  ) 
r'r-  -  G  N  J  l*  ^  G  :  +  r,  \  •  7  *  ? 

T : c 
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jPA«.c  >  n 

oc.s  *  no 

19  IFCOG.GE.DGSJfO  TO  9 

DCS  *  OG 
XPS  -  XP 
Y P  S  *  TP 


XPP5  =x°p 
yppi.yoc 
IP  5  *  2  P 

jpflts  ZPB 

9  IF{ID.Nf.01W9!TP<6»l,'lX(!) 

10  FOR1*1  &  T  (  1  h  9F 1 1 .  ?  »f  1  3  •  4  1 

11  Zp  *  Z  P  *  PFlY* 
JFtZP.Lf.^'GO  TP  1? 
xp  s  xP  +  n p L T  t 

1FJXP.LT. XP^AXI'O  TO  8 

J3  XPP  =  XPP  ♦  OFlTA 


y ( n .z 1 1 1 .xp.tp.zp.xpp  .vpp  ,zpp 


1F(X0P.L5:.0.  >G0  TO  5  „ 

*RtTE(6*10)Xt!).vil)*Z<n»xPS«''’P'» 

14  CONTINUE 

15  C-0  TO  2 
ENO 


zPs,xpps.TpPe*«zppr 


DOG 


4b 


IIRFTf  X-2-A  FUL ! ST 

DIMENSION  X  (  50  ) *  Y ( 50 ) »Z ! 50 ) .TITLE! 13) .VACV I 5 ) . V8C V ( 5  I »VACU<5>  » 

1  VPCU( 5 1 .RACV(5>  .RACUI5I . RBCVI5) iRBCUIS ) 

NAMELIST/ I NPUT/  A.B.C.D.E.f.R.T.DELTA.fpsP  »EPSPp»DEACV  .OEBCV. 
1DEACU.DEBCU.VACV ,VBCV.  VACU.VPCL  ACV.RRCv  .RACU  «RBCU  *  X . Y ,2 

1  FORMAT ( 1H139X53HA  EROjET-GENERal  CCRPORATI 
10  N/1H039X53HC  0  M  P  U  T  I  N  G  SCIENCES  DIVISION 
2/1 H050X3 1 H3  ACRAMENTC  PLAN  T  //IHC57X17Hr‘  R  0  G  R  A  V 

3  -C/1HC19X93HG  FNERAl  TwC-°hASE  EOuIlIBRI 

4  U  M  IN  TERNARY  SYSTE  M/16M0C  S  T  0  M  E  R  97  x 

519HP  ROGRAMME  R/12H0Y.  A.  CH ANG  1  09 X 1 1  HR .h.  M4RlER//) 

RE  AD  I  5  *  3  )  TIT'. F. ID 

3  FORMAT  11346,17  ) 

2  REAP  (5. INPUT) 

WRITF  (6.1) 

WRITE(fe»4)TITLE.IP 

4  FORMAT ( 1H026X13A6.25X  I  2  ) 

WRITE  (6. INPUT) 

WRITEI6.1 ) 

WR I TE ( 6  *  1 6 ) TITLE 

16  FORMAT  (  1H02  6X 1  3  A  6  /  /  7xlHXl2XlnYl2XlH2nx2HXPllX2HYPllX2HZPllX3HXPP 
1 1 l X3HYPP1 1X3H2PP  9X7H0ELTA  0/ ) 

flag  =id/io 

BT  =  R*T 
DO  14  I  =  1.50 

I F ( x ( I  )  .FQ.O. . AMP.Y ( ! 1 .FG.O. .ANO.Z ( I ) .FQ.n. )G0  TO  15 
I  PASS  =  1 
XPP  =  A 

00  17  IXPP=l.lOO0 
I  F  (  X  PP  — X  (  I  )  1  6  ♦  1  ■*  »  7 

6  XP  =  Xi  I  )  4-  DELTA 

X  PM  A  X  =  F 

GO  TO  8 

7  XP  =  F 
XPVAX  =  X  t I) 

8  DO  lb  IXP  =  1.10P'' 

ZP  =  C 

pp  ??  I ?P  = 1  .  1  OOF 
YP=) ,-XP-ZP 
! F ( XP . EO. XPP ) GO  TO  11 
GNU  I  =  <X(  D-XPP  t /(XP-XPP) 

IF(GNUl.3E.1..0R.GNUl.LT.O.)GC  TO  11 
GNU2  =  1 •-GNU! 

Y  PP  =<Y( I)-oNUl*YP)/GNU? 

Z  PP  =1,-XPP-YPP 

I F (ZPP.GF.l ..OR.ZPP.LT  .0 . ) GO  TC  11 

CALL  DLGC ( DGACV ,R ♦ T .ZPP . VACV . R AC  V . DE AC  V .F L AG . I 1  1  ) 

CALL  DLGC ( DGRC  V  »R ♦ T .ZPP.VACV.RSCV.DEBCV.FLAG.Ill) 

U 1  =  XPP/  (  1  «-.?  PP  t 
IJ2  =  Y PP/  (  1  .-ZPP  ) 

CALL  XLNXIU1.FU1  .Ill  ) 

CALL  XLNX CU2 .Fu2 .51 1 ) 

DG2=  (  (  Fjl  +FU2  )  *PTfUl  *U2*EPSPP  +  L'2*OG3CV  +  Ul*OGArv  )*U  .-ZPP  ) 

CALL  D'_GC  (  OGACU.R.T  «  Z  n  .VACU.RACU  .DE  ACU  .  FL  AG  .  5  1  1  ! 

C  A.  L  OLGCOGBCU.P.T  .Zp.VRC.  DEBCu.Fl  a-G.II  1  ) 

Ul=  XP/11.-ZP1 
i.z=  '  :  .-op  i 

call  x i  *.x  c  ;i  1  . r ■  j i  ."ii  ) 

CALL  C.2,F-.:;  .Ill! 
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19 


DG1  =<  tFul+FV?  l*RT*Ul»V2«EPr-P*v!? 
DO  •  GNul*DGl  ♦  0"yU2  •  r'G? 

1  P(  iPAX.^.m.'-'  ir.p  rr>  tq 
IPm  *  o 
00  E  i  D"'- 

1  F  (  OG •  G t  •  C  G  9  GO  T  c ,  "* 

JGS  *  Of, 

XfS  =  XF 
YPS  «  Y  P 


:u+Gl  *GGxr'.n  •  1 1 .  — zr* » 


XPPS  *APP 
VPPXcYPP 
?P<;  s  /P 
2PP<;=  Zop 

9  I  F  [  10.  nE  .01«.'*ITr,6.i-|Y  ,  ;  j  :  :  .,YPtvot;  =  tXsP.vPr 

10  FORMAT  (  1*  9P  1  1.-'  «c  1''.- ) 

1 1  ZP  =  IP  -  DEL  *  * 

If  tZP.GT.'.lC-O  TO  12 
21  CUNT  !  VUF. 
wR  1  TE  t 6.2  3) 

21  FORMAT  (  25-0  T->E  lGC"  .•Ar  *2  x  0 E  r -E  0  ) 

GO  TO  5 

12  *?  SXP+DEL  T  A 

1  FIXP.Gf  .XP'’AX  )'C.  rr'  IX 


P».D 


18  CONTINUE 

WPlTf I  6.20  I 

2*  FORMAT  (  2  5*-pth*  «[•  LOOP  •:•*<  : 


r.O  TO  5 

1 5  XPP  *  XPP  ♦  DELTA 
ir ( XPP.GT.r IGO  TO  5 
1 7  r ONT 1NDE 
WRITE (6.21  ) 

21  FORMAT  1 26H0THE  XPP  LOG'-’  A  0  E  x  0  L  F^E  D ) 

5  W«  I  T  E  (  6 . 1 C  )  X  (  )  .  >  t  I  i  .  2  I  ;  )  .  X  p  •  .  Y P  f. ,  o  , 

14  CONTI Not 

15  GC  TO  2 
F  NO 


,  Y  P  P  s 


4b 


r 


i  v  ‘  F  i  L  I  S  T 

P  I  vrf.jc;  I0N|  TULrini 

HMMON  COGACV  (  3  )  ♦  C  DGBC  V  (  3  )  ♦CDGACW(  3  )  ♦CDGBCWf  3  ) 

NAVr  L  I  ST  /  I  NPtlT  /  Pi  ,13  2  .DELTA  .X  ,  Y  ,Z  .XP  ,YP  ,ZP  tXPP.YPP  .ZPP.R  ,T  fM, 
KDGACV  ♦  C'PGBCv  .  COGAC  w  ♦  CDGBCW  *  EP  RP  *EP  SPP  .EPSPPP 

1  iv-KMAT  t  lHi  j'JXP  "MiA  (.  RO  JET-GENERAL  CORPORATI 
*v 1H039X5  3MC  0  V  P  U  T  I  N  G  SCIENCES  DIVISION 
v  /  IhCbGXMHj  A  c  R  A  V  t  N  T  0  PLAN  T / / 1 H05 7X 1 7HP  h  0  G  R  A  M 
3  K  /  1 H  0  ’  X  9  3  H  G  E  N  E  R  A  L  TWO-PHASE  EQUILIBRI 

^  B  o  l  r  n  T  E  R  N  A  R  Y  5  V  S  I  E  M/16HOC  U  S  T  0  M  E  R  9)X 

blVMP  R  0  G  R  A  '•  M  L  R/12HOY.  A.  CHANG109X11HR.H.  MARLER//) 

PEADIS.3TT ITLE .ID 
3  FORMAT (  l 3 A 6  *  I  3  ) 


3  R  F  AD  (  S  ,  INPUT  ) 

WRITE  (6.1) 

WR ITE (6  *A ) T I TLE  »  10 
A  FORMAT (  1H036X  1  3A6, 25X12  ) 
WRITF  (6. INPUT) 

WR I TE 16 . 1 ) 

W  R I T  E ( 6  *  5 ) TITLE 


w  r  =  r  »  t 

CALL  DGACVS1DGACV.2PP 
CALL  DGBCVS1DGRCV.2PP 
Ml=XP/( 1.-2P) 

1  j  ?  =  Y  P / ( } .  -  7  P  ) 

CALL  XLNXdJl  ,  FIJI  ,  Til  ) 
call  xlnx UJ? .Fll? ,«  1  1  ) 
PG  ]  =  (  F  U  1  +  T  U  3  )  *  R  T 

i ;  1 =XPP/ ( 1 . -ZPP ) 

U  3  =  Y  P  P / ( 1 . -ZPP ) 

CALL  XLNX ( U1 »FUI ♦ SI 1 ) 
CALL  XLNXCJP.FUZ  .SI  1  ) 
DG2  =  ( ( F J1 +  FU? ) *RT 

1  1  I  PASS -1 


♦  111  ) 
♦  SI  1  ) 


+  ERSP*U1*1J2 


+EPSPP»U1 *02 


+  DGBCV*U2  +  DGACV*U1)*(1 


-ZPP  ) 


HYP  -  YP-Y 


PYPP= YPP-Y 
P7P  =  Z  P-7 
~7  D  n -  7  D  n  —  7 

A p  F  A  3  =r')YPP*07P  -  0VP*3ZPP 

Z  P  P  d  =  p  i 
7  Y  P  P  P  =  0  . 

Y  D  P  D  M  X  =  ]  ,-ZPPP 
n  2Dop=7.  PPP-7. 

P  f: YPPP  =  yppp-y 

A  < E A  1  =DYP«JZPPP-DYPPP*DZP 
A.REA3  =DYPPP»C2PP  -  DYPP*DZPPP 

I E ( . NOT . ( AREA  1 .LT . 0. . AND.AREA2.  L  T.O. . AND. AREA3 .LT.O.  .OR. AREA1 .GT . 
1  r. .  .  an:  .ARE  A2.GT  .0.  .  AND.AREA3.GT.0.  )  )G0  TO  11 

X  P  P  P  = 1  .-YPPP-7PPP 

CALL  HGACWS ( 3GACW ,ZPPP  .Til) 

CALL  OGRCWSIDGP-W.ZPPP.J 1 1  ) 

I)  I  =XPPP/  ( 1  .-7PPP  ) 

1  ..'3  =  YPPP/  {  1  .  -ZPP p  ) 

CALL  XLNX ( O]  .Ftn  ,  S  1  1  ) 

CALL  XLNX ( U? .PJ2  .  U  1  ! 

303  =  (  (Ful+Fj2)«RT+t"P5PPP*'Jl«U2  +  DGBCW*U2  +  DGAC W*U 1  )  *  (  1 . -Z PPP  ) 
l)  =  (rPp  -  Y  p  )  *  ( Z  P  P  P  -  ZP)  -  (ZPP  -  ZP)  *  (YPPP-  YP) 
i ) (D.EG.n. ) GO  TO  1 1 

G  ‘it;  1  ^  ((  YPP-Y  )»(  Z°PP-Z  i  -  (  ZPP-Z  )*  (YPPP-Y  ))  /D 
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IF<6NU1.LE.0..0R.GNU1.GT.1.)60  TO  11 

GNU2  ■  I IY-YP)*IZPPP-ZP>  -  (Z-ZP)*IYPPP  YP)J/D 

1F(6NU2»LE»0«*0R«GNU2 •GT«1» >60  TO  11 

GNU3  *1.-GNU1-GNU2 

06  •  GNUl*  0G1  +  GNU2  *DG2  GNU3&D63 
IFUPAS$.EQ.O>60  TO  9 
I  PASS  ■  0 
DOS  -  06 

9  I F (DG»Gc»DGS)GO  TO  10 
OGS  •  OG 
XPPPS  »  XPPP 
YPPPS  *  YPPP 

to  imL!U"”mE(6.1MX,Y,Z.XP.YP.ZP.XPP.YPP,2l>l>.Xf>PP.YPPI>.ZPI>(>.C>6 

11  YPPP  «  YPPP  ♦  DELTA 
IF(YPPP.LE.YPPPMX>GO  TO  8 
ZPPP  •  ZPPP  ♦  DELTA 
IF(ZPPP.LE.B2>  GO  TO  7 

WRITE«6»14>X.Y»Z»XP  »YP  *ZP  »XPP  »YPP  *ZPP  * 

1 XPPPS » YPPPS *ZPPPS  »DGS 

M  *  M  -1 

IF (H»LT *0)60  TO  2 
IF(DGS.GE.0.»GO  TO  12 
B1  -AMAXUP1.0GS-DELTA5 
B2  «AMIN1CG2*DGS+DELTA> 

12  DELTA  ■  DELTA/ 10 • 

GO  TO  13 

14  FORMAT ( 1H  12F10.7  »F10»3 ) 

5  FORMAT C1H026X  _  xP  Y 

1  13A6//130H  X  Y  TDD  Z  yPDp  Yppp 

ip  zp  XPP  YPP  ZPP  XPPP  YPPP 

2  ZPPP  DELTA  G/> 

END 


YPPP 


SO 


I 


n 


SIBFTC  K-3-A  FULlST 

DIMENSION  TITLE (13) ,PACV<5) .RBCV<5) ,RACW!5 ) .RBCW!5  > .VACV15) t 
1  VBCV ! 5 ) » VACW ! 5 ) *  VBCW ( 5 ) 

NAMELI5T/INPUT/B1,62*DELTA,X,Y ,Z ,XP »YP ,ZP »XPP. YPP »ZPP .R » T .M  .EPSP ♦ 

1  EPSPP  .EPSPPP  .DEACV .DERCV .DEACW .DEBCW. VACV  .VBCV .VACW  »VBCV,' » 

2  RACV.RBCV.paCW.RBCW.DGAOJ.DGRCU 

1  FORMAT! 1HX39X53HA  EROJET-GENERAL  CORPORATI 
10  N/1H039X53HC  0MPUTIN6  SCIENCES  DIVISION 
2 '1H050X31HS  ACRAVENTO  °  L  A  N  T//1H057X17HP  R  0  G  R  A  M 


3  K/1H019X93HG  ENERAL  TWO-PHASE  EQUILlBRI 

4  U  M  IN  TERNARY  SYSTE  M/16HQC  U  S  T  0  M  E  R  97X 
519HP  ROGRAMME  R/12H0Y.  A*  CHANG109X11HR.H.  MARLER//) 

RFAD!5.3)T ITLE. ID 
5  FORMAT ( 13A6* 12 ) 

>  READ  ( 5  » INPUT ) 

WRITE  ( 6 1 1 1 
WRITEI6.4J TITLE. 10 
4  FOR“ATt 1H026X13A6. 25X12) 

WRITE  (6.  INPUT) 

WRITF16.1 ) 

WR ITE(6»5)TITLC 
rL AG  =  ID/10 
R  T  s°*T 

CALL  DLGCIDGACV.R.T. ZPP.VACV.RACV.DEACV.FLAG.Sll) 

CALI.  DLGCfDGBCV.R.T.ZPP.VBCV, RBCV.DEACV.FLAG.Sll) 

■J1«XP/(1.-ZP) 

U2=YP/<  ,-ZP) 

C*LL  XLNXIUl.FIJI  .Sll) 

CALL  XLNXt')2.FtJ2*«ll  ) 

DG1  =  !cUl+FU2)*RT  +  epcp*tji*i i2+DOACU*Xp-*-DG3CU*vP 
U1  *X°P/ ( 1 . -7Pp) 

U2=YpP/  ( 1  • -ZnR 1 
CALL  XLNXtUl. FIJI. Sll) 

CALL  XLNX(,J?.r,J2.Sll  ) 

DG2  =  (  (  FU1+CU2  )*RT  +  EpSPP*Ul*tJ2  ♦  DC,BCV»U2  +  DG  *CV*U1  )  *  (  1 .  -ZPP  1 

3  IPASS=1 
DYP  =  YP-Y 
pyoosYPP-Y 
DZP  =  ZD-Z 

ozpp=zpp-z 

AREA?  =DYP°*DZp  -  DYP*DZPr> 

ZPPD=  OI 
7  yppD  =  o. 

Y?ppmx=1,-ZPP° 

^ZPPP-ZPPP-Z 
9  DYPno=vppp_y 

AREA1  *DYP»DZPc,r'-nYPP®*DZn 
AREA3  =DYPPP*DZ°P  -  CYpr«*DZpPD 

I F ( .NOT . t AREA l.LT .0«. ANO. ARE \2.'-T.0». AND. ^RE"3.LT.0..0R.AREA1.CT. 
10..AND.AREA2.CT.0..AND.AREA3.GT.C. ) )GD  TO  11 
XPPP=1 .-yopP-ZPpP 

CALL  DlGC ( DGACW.R .T .ZP°P .VACW . RACW . DE ACW . F LAG .51 1 ) 

C "LL  DLGC! DG’CW.o.T .ZPPP  .  VRCW  ."BCVt .  DEPCW  »FL  AG » «  1 1 1 
•jl=Xp=>°/(l.-ZPOC') 

!l?=voop/  ( 1  .-zonn  ) 

CALL  XLNXIU1.FU1 .Sll) 

CALL  XLNXri2.F  I2.su  ) 

DG3  =  11FU1+FU2  *RT+rPr.PPP»'Jl*U2  +  DGBCW*'J?  +  DGACW*’J1 )  *  ( 1  .-ZPPP ) 
0  =  ( YP°  -  Y°)*(ZPPP-  ZP>  -  !  ZPP  -  ZP)  *  ;yppt>_  vP) 


c-pco*tji*u2+D0ACU»Xp-fDGPCU#vP 


DGACV*U1 )*! l.-ZPP) 
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ifco.eo.o*  >oo  to  11 

GNU1  ■  ( <VPP-Y)*(ZPPP-Z|  -  <ZPP-2)*(YPPP-Y)|/D 

IF(GNUl.LE.O**OR.GNUl.GT.l.)GO  TO  11 

GNU2  ■  { (Y-YP|*{2PPP-ZP»  -  ( 2-2P )* < YPPP-YP ) ) /D 

IFIGNU2.LE.O..OR.GNU2.GT.1.1GO  TO  11 

GNU3  ■1.-GNU1-GNU2 

OG  •  GNUl*  DG1  +  GNU2  *DG2  +  GNU3*D03 
IFIIPASS.EG.OlGO  TO  9 
IPASS  ■  0 
DGS  •  DG 

9  IF (DG.GT.DG5IG0  TO  10 
DGS  *  DG 
XPPPS  *  XPPP 
YPPPS  r.  YPPP 
2PPP5  ■  ZPPP 

10  I F ( ID.NE.O) WRITE  I  & • 14  I X  »Y»Z»XP» YP  *ZP*XPP»YPP» ZPP» XPPP  * YPPP  #2PPP»DG 

11  YPPP  «  YPPP  +  DELTA 
!F<YPPP.LE.YPPPMX)G0  TO  8 
2PPP  *  ZPPP  DELTA 

IF ( ZPPP. LE. 92 1  GO  TO  7 

WR1TE(G«14)X»Y*Z*XP  »YP  »ZP  .XPP  .YPP  ,ZPP  t 
1XPPPS. YPPPS .ZPPPS.RGS 
M  ■  M  -1 

IFIM.LT.OIGO  TO  2 
IF(DGS.GE.O.)GO  TO  12 
B1  -AMAX1 <  Bl.DGS-DELTA J 
B2  »AMINUB2.0GS+DELTA) 

12  DELTA  •  DELTA/10* 

GO  TO  11 

1*  FORMAT ( 1H  12F10*7  »F10.3) 

5  FORMAT (1H026X 


1 

13A6//130H 

X 

Y 

Z 

XP  Y 

IP 

2  ZP®P 

END 

ZP  XPP 

delta  g/) 

YPP 

ZPP 

XPPP 

YPPP 

52 


£ 

fi 

I! 

0 

0 

fl 


94  8  93  77  77496405 
'  7  9  5  Z  *)P0 
7  *G  074074*  74 •449* 


$  I BLDR  XLNXF 
STEXT  XLNXF 

*N  < <XS7v*7'*7-5  00  80  9  9 

*N9( ( * *XV*7 ( *7— 5  5  •  74+03*  W 

*N8V*'  *7V*X 1*7-74  04*74774  •  *5 

*59(2'*/  •  *5  •  • 

SCGICT  XLNXF 

*/  ()(  0*  l'*)P  0  0  9* ) P  '*>P-6  P*J- 

SDKEND  XLNXF 

SIBFTC  CUBICF  NOLlST 

SUBROUTINE  CUBIC(COEF.X»Y.*) 

DIMENSION  COEF ( 3 ) 

IF(X.EQ.l. JRETURN  1 

Y  *  ( (X*COEF(3)+COEF<2> )*X+C0EF(  1) ) /(l.-X) 

RETURN 

END 


12/17/65  XLNXOOOO 
12/17/65  XLNX0001 
Z4  8  93  *6  *-5  *  XLNX0002 
• ) 4*  G47* 4*  04* 74 *  XLNX0003 
*•  7*  *•  7447*  • XLNX0004 

XLNX0005 
12/17/65  XLNX0006 
XLNX0007 
12/17/65  XLNX0008 


G 

I 

I 

II 

[ 

I 

£ 

L1 

| 

I 
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SIBFTC  DGACVF  FULIST 

SUBROUT  INF  DGACVS!D6ACV#ZPP  **> 

COMMON  CDGACVO  )  .CDGBCV  O  >  .CDGACWt  3 )  .CDGBCW 1 3 ) 

CALL  CUBIC!CDGACV,ZPP.DGACV»S10> 

RETURN 
10  RETURN  1 
END 

S1BFTC  DGBCVF  FULIST 

SUBROUTINE  DGBCV5{DGBCV*ZPP  »*> 

COMMON  CDGACV ( 3  >  »t DGBCVt  3 ) .CDGAC W< 3 ) • CDGBCW ( 3 ) 

CALL  CUBICICDGBCV.ZPP.DGBCV.SIOJ 
RETURN 
10  RETURN  1 
END 

♦I8FTC  OGACWF  FULIST 

SUBROUTINE  DGACWSt DGACWtZPPP.* ) 

COMMON  CDGACVO). CDGBCVO)  .CDGACWO)  .CDGBCWO) 

CALL  CUB I C ( CDGACW *ZPPP  *DGACW  »  S 10 ) 

RETURN 
10  RETURN  1 
END 

SIBFTC  DGBCWF  FULIST 

SUBROUTINE  DGBCWS! DGBCW.ZPPP.*) 

COMMON  CDGACVO  )  .CDGBCVO)  .CDGACWO )  *  CDGBCW  (  3) 

CALL  CUBIC! CDGBCW »ZPPP .DGBCW * S10 ) 

RETURN 
10  RETURN  1 
END 
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SIBFTC  DELGC  PULIST 

SUBROUTINE  DLGC  ( DGC  .R  »  T » ARG  W.  ARE  .DELE  .FLAG  .* ) 

DIMENSION  V<  5 } .ARE ( 5  ) 

GC  •  ARE ( 5  )  *T  ♦  ARE ( 4 ) 

DSFC  ■ ( ALOGIO ( T ) *ARE ( 2 ) +ARE { 3 1 ) *T  +  ARE ( 1  1 
CALL  DGXSl ARG.DELE .R .T .DELGXS .FLAG .S2 ) 

DGC«( ( I <ARG*V(5 >+V<4> )*ARG+V< 3  >  )*ARG+V<2) )*ARG+V< 1 )+( l.-2.*ARG)* 
1GC-M  1 .-ARG ) *DGFC+DELGXS ) / 1 1 .-ARG ) 

IF (FLAG.GT.O. 1WRITE ( 6*1 1GC .DGFC.DGC 

1  FORMAT (  6H  GC  *  E15.8.5X  7HDGFC  «  E15.8.5X6HDGC  «  E15.8J 
RETURN 

2  RETU9N1 
END 

SIBFTC  DlGXS  FULIST 

SUBROUT  I NE  DGXS ( X .DELE .R.T .DELGXS. FLAG . * > 

RT«R*T 

I F (DELE) 10.5*1 
1  W-DELE/RT 
EXPW  •  EXP ( W ) 

EXPNW  -l./EXPW 

CALL  XSUBABI X.EXPNW.XA.XB.FLAG.S10) 

CALL  PRELIMtXA.XB.EXPNW. EXPW. SMALLR.S. FLAG) 

CA LL  DELGCX ( DELE .RT . XA . XB. SMALLR.S . DLGCXS. FLAG ) 

CALL  DELGMXIXtXA.XB.RT .DELE. DLGCXS. DGMEXS. FLAG) 

CONJXA-l.-XA 
C0NJXB-1.-XB 
CALL  XLNX(XA.FXA.SIO) 

CALL  XLNX(CONJXA.FCNJXA.SIO) 

CALL  XLNX(XB»FXB»S10) 

CALL  XLNX(CONJXBtFCNJXB.SIO) 

DELGXS- 1  (FXA+FXB+FCN JXA  ’•FCNJXB  )*RT  +XB*DELE ) *U  .-X  )  /2  . 

GO  TO  9 
5  CONJX-l.-X 

CALL  XLNX  < X.FX.S10 ) 

CALL  XLNX  (CONJX.FCONJX.SIO) 

UX*1.-2.*X 

CALL  XLNX(UX.FUX.SIO) 

DELGXS  * ( FUX-FCONJX+  FX)*RT 
9  I F (FLAG. GT. 0. ) WRITE l 6. 11) DELGXS 
RETURN 

10  RETURN  1 

11  FORMAT ( 10H  DELGXS  -  E15.8) 

END 
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•EXECUTE  IR  JOB 

SIBJOB  A4B052  MAP 
SIBFTC  TAB  FULI ST 

DIMENSION  DELTAE 110 ) »TEE< 10 ) 

1  READ! 5»100I IDE* IT* (DELTAEI I J » 1 «1 *  IDE) *  < TEE ( I ) »I » I T ) »XF . XL .DX.R . 
1FLA6 

WRITE (6 *500 1 <DELTAE( 1 )»!*1»IDE>» (TEE(I)»I«1*IT) .XF.XL.DX.R.FLAG 

DO  50  l-l *  IT 

T-TEEI1) 

RT«T*R 

DO  49  J»1 » IDE 
DELE-OELTAE(J) 

WRITE  Ift.2001T.DELE 

W-OELE/RT 

EXPW  ■  EXP (W) 

EXPNW  ■  EXP(-W) 

X«XF 

lFtFLAG.GT.O.)  WRI TE l 6*300) RT *W> EXPW. EXPNW *X 

10  CALL  XSUBAB(X.EXPNW.XA.XB.FLAG.Sll) 

CALL  PREL I M ( X A  *XB .EXPNW  *  EXPW  * SMALLR  *S • FLAG ) 

CALL  DELGCX { DELE *RT »XA • XB. SMALLR *S * DLGCXS .FLAG ) 

CALL  DELGMX IX. XA.XB.R., DELE. DLGCXS. DGMEXS. FLAG) 

DELGXS*  ( ( X A*ALOG  ( XA ) «-XB*ALOG l XB  > ♦ 1 1 . -XA )  *ALOG 1 1 .-XA  >  +  ( 1 . -XB ) *ALOG ( 

1 1 .-xb  n *rt*xb*dfle  >  *  n .-x ) n . 

WR I TE( 6. *00 )X. DLGCXS. DGMEXS. DELGXS 
I F ( FLAG.EO. 1. )FLAG«0. 

11  XmX+DX 

IFIX.LE.XLIGO  TO  10 
*9  CONTINUE 
50  CONTINUE 

IOC  FORMAT 12I2/(5E13.7  )) 

200  FORMAT C1H139XAHT  *  E15.8 » 10X10HDELT A  E  »  E15.8/1H038X1HX14X12HDELT 
1A  G  C  XS  7X13MDELTA  G  ME  XS  9X10HDELTA  G  XS/ ) 

3 Of  FORMAT  1 1H026X9E16.8 ) 

*00  FORMAT ( 2TXAE20.8 ) 

500  FORMAT ( 1H1/I 26X5E1G.8  )) 

GO  TO  l 


* IBFTC  XSBAB  FULIST  A 

SUBROUTINE  XSUBAB!X.EXPNW.XA.XB.FLAG.*>  B 

ASTAR*1.-EXPNW  C 

CONJX*l«— X  D 

8STAR* ( ( 1.+X)*£XPNW-3.*X+1. J/CONJX 
CSTAR«-2.«X/CONJX*EXPNW 

XB1—BSTAR/2. /ASTAR  G 

XB2*SQRT !BSYAR*BSTAR-4.*ASTAR*CSTAR ) /2 • /ASTAR  H 

XB-XB1+XB2  I 

IF!  (XB.GT.O.I.AND.IXB.LT.l.nGO  TO  10  J 

XB-XB1-XB2  K 

IF! (XB.GT.O. > .AND. (XB.LT.l. >  )GO  TO  10  L 

WRITE  (6.100JXB. ASTAR. BSTAR.CSTAR  M 

100  FORMAT  1 6H0XB  «=  E15.8.38H  .  IS  NOT  IN  THE  RANGE  10  LT  XB  LT  1)./  N 

1  9H0ASTAR  -  E15.8.10X  fiHBSTAR  *  E15.8.10X8HCSTAR  «  E15.8/  0 

216HOGO  TO  NEXT  CASE)  P 

9  RETURN  1  Q 

10  XA*«2./C0NJX*X-XB  R 

IFUXA.GT.O.  )  .AND.IXA.LT.  1.  ))G0  TO  12 

WRITE  !6.200>XA. ASTAR. BSTAR.XB  T 

200  FORMAT  16H0XA  *  E15.8.38H  .  IS  NOT  IN  THE  RANGE  !0  LT  XA  LT  H./  U 

1  9HOASTAR  =  E15.8.10X  8HBSTAR  =  E15.8.10X  8HCSTAR  =  E15.8.10X  A 

25HXB  *  E15.8/16HOGO  TO  NEXT  CASE)  B 

GO  TO  9  C 

12  IF (FLAG.LE.O. )GO  TO  1*  D 

WRITE  I  6.300 JXA.XB .ASTAR .BSTAR.CSTAR  E 

300  FORMAT!  21H0D  UMP  XSUBA  B/6H0XA  =  E15.8.5X5HXB  *  E15.8.5X  F 

18HASTAR  =  E15.8.5X8HBSTAR  =  E15.8.5X.8HC5TAR  *  E15.8/)  G 

1 A  RETURN  H 

END  I 

SIBFTC  PRLIM  FULIST  A 

SUBROUTINE  PRELIM!  XA.XB.EXPNW.EXPW.SMALLR.S.FLAG)  B 

CONJXA-1 .-XA  C 

CONJXB=l«— XB  D 

U=XA /CON JX A*EXPNW  E 

V=XB/CONJXB*EXPW  F 

RNUM-! !l.+U)*CONJXA)**2  G 

SMALLR=RNUM/!RNUM+EXPNW>  H 

SNUM= ! 1 1 .+V ) *CONJXB ) **2  I 

S*SNUM/!SNUM+EXPW)  J 

IF«FLAG.LE.O. )GO  TO  1  K 

WRITE  !6.100)U»V»RNUM,SMALLR.SNUM,S  L 

100  FORMAT ( 21HOD  UMP  P  R  E  L  I  M/5HOU  =  E15.7.5H  V  =  E15.7.8H  RNUM  M 

1=  E15.7.  9HSMALLR  =  E15.7.8H  SNUM  =  E15.7.5H  S  =  E15.7/J  N 

1  RETURN  0 

FND  P 
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SIBFTC  DGCX5  FULIST 

SUBROUTINE  DELGCX  ( DELE .RT ,XA  <X8 * SMALL R » 5  .DLGCXS.FLAG  > 
OLGCX5»(ALOG(XA/U#-XA>  l*SMALLR  *  A  LOG  l  X  8/  { 1.-X8  J  }*S) #RT  ♦  S*DELE 
IFIFLAG.GT.O. J MR I TE ( 6 . 100) DLGCXS 
RETURN 

100  FORMAT ( 21H0D  UMP  D  L  G  C  X  S/13H0  DEL  GCX5  «  E15.7/I 
END 

SIBFTC  DGMX  FUL 1ST 

SUBROUTINE  DELGMX  CX.XA.X8 iRT  .DELE .DLGCXS .DGMEXS .FLAG) 

CONJXA-l.-XA 

C0NJXB«1.-XB 

Z AMX« ALOG I  CON JXA  J  *CON JXA  ♦  ALOGCXA)*XA 

ZBMX«ALOGICONJXBI*CONJXB  ♦  ALOG( XB ) *X8 

DGMEXS- ( <  ZAMX+ZBMX ) *RT+XB*DELE > / 2. -X/ I 1 .-X  ) *DLGCXS 

IF (FLAG. GT .0* IWRITE  <6*100)ZAMX.ZBMX,tGMEXS 

RETURN 

100  FORMAT (21H00  UMP  DELGM  X/8H0ZAMX  *  E15.7.10X7HZBMX  ■  E15.7, 
110X9HDGMEXS  »  El 5. 7/) 

END 

SIBFTC  OGMS  FULIST 

SUBROUT I NE  DOMES ( X • RT ♦ ALPHA .GME * DGMESG . GC tDELGF .FLAG ) 

IF(X-. 5)10. 20.30 

10  DGMESG  »ALOG(  Cl.-2.*X)/(l.-Xn*RT  ♦  2,*DELGF  ♦  GC 
GO  TO  40 

20  DGMESG  •  ALOG( (l.-2.*ALPHA»/ALPHA)*RT  -  GME 
GO  TO  40 

30  DGMESG  »  ALOGC  U.-X)/(2.*X-1.H*RT  -  GME 
40  IF(FLAG.GT.O. )WRITE(6.100J  DGMESG 
100  FORMAT ( I 9H0D  UMP  DOME  S/IOHODGMESG  =  E15.B/I 
RETURN 
END 

SIBFTC  DGCS  FULIST 

SUBROUTINE  DGC5IGIX.GC.GME .ALPHA, DELGF.RT .DLGCSG. FLAG! 

IFfX-.5llC. 20.30 

10  DLGCSG  «  AL0G(X/(1  ,-2.*X»)*RT  -GC 
GO  TO  40 

20  DLGCSG  «  ALOG ( ( 1 .-2 ,* ALPHA ) /ALPHA  J  *RT -GC 
GO  TO  40 

30  DLGCSG  «  ALOGC (2.»X«1.)/X)»RT  +  2.*DELGF  ♦  GME 


40  IFIFLAG.GT.O. J WRITE (6. 100) DLGCSG  1 

100  FORMAT(21HOD  UMP  D  G  C  S  I  G/10H0DLGCSG  *  E15.7/J  j 

RETURN  K 

END  L 
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APPENDIX  II. 


Computer  program  for  calculating  the  excess  free  energies 
of  the  M«2C  phases  in  terms  of  2  +  w  gramatom  of  alloy.  For 
theoretical  background  of  this  problem,  refer  to  AFML-TR-65-2, 
Part  I,  Volume  I  {1965}.  ........  (Bill  Reuse). 


59 
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*6X£CUTE  JBJ08 

S I e JOB  42051  MAP 

SIBFTC  42051  FUL! 5T .DECK < REF 

C  PROGRAM  C  (42051)  ME2CY  PHASE  CALCULATION  -  EXCESS  QUANTITIES 
DIMENSION  Y ( 2500 ) *  XCA(2500)  .  XCBI2500J.  ZA(2500)*  ZB(2500). 

1  2(2500).  MXSI25005.  GXS(2500) 

ID  READ  (5*20)  MODEL*  YMI N»  YMAX.  DY.  WMIN*  WMAX,  D W 
.  FORMAT  ( II .8X.6E10.0) 

WRITE  (6*30)  M0DEL»YM1N*YMAX*DY*WMIN»WMAX.DW 
30  FORMAT  (lHl//////<'/////////61X^9HPR0GRAM  C/ /62X  »7H  ( 42051 )  /  /  /44X  . 

1  44HME2CY  PHASE  CALCULATIONS  -  EXCESS  QUANT  I TI ES/ //62 X .6HMODEL  » 

2  I1////3X.6HYMIN  «  F11.5  ,4X,  6HYMAX  «  Ell. 5  *4X.  6H  OY  *  F11.5* 

3  4X.  CHWMIN  *  F11.5.4X*  6HWMAX  *  F11.5  *4X.  6H  DW  »  F11.5  ) 

W  *  wmin 

35  LINCT  «  50 
37  DO  79  1*1  *2500 

Y  (  I  )  *  0. 

XCA( I )  *  0. 

XCB( I )  *  0. 

ZA (  I )  -  0. 

ZB ( 1 1  «  0. 

2(1)  *  0. 

HXS(I)  ■  0. 

39  GXS(I)  *  0. 

40  I  *  1 

Y ( 1 )  «  YMIN 

EMW  *  1./  FXP(W) 

50  GO  TO  (100*200) .MODEL 

100  AA  *  2.*  (l.-EMW) 

AAAA  *  2.  *  AA 

105  B  -  1.  -  Y(I)  +  EMW*( l.+YI I ) ) 

C  *  Y (I)  *  EMW 
OU AO  *  P**2  +  aaaa*C 
IF  (QUAD)  107*110.110 
107  24(1)  =  l. 

ZP(I)  *  QUAD 
GO  TO 

no  quad  *  sort ( quad ) 

XCBC I )  *  ( -P+QUAD ) / AA 
IF  (XCB( I )-l* ) 120*130*130 
120  7  F  (XCB(I) ) 130.130*160 
130  ZB(  I  )  *  XCB (  I  ) 

\  XCBU)  *  (-B-QUADJ/AA 

IF  (XC“(I))  150.155*140 
140  IF  |XC°.(  H-l.)16'‘. 150*150 
1  cr>  2  (  t  >  =  X"  ( i ) 

XC°< 1 )  =0. 

155.ZM1)  =  3. 

GO  jn 

160  XCA  (I)  *  Y<  n  -  XCB (  I ) 

ZA(I)  =  XC  A  (  I  )  *  ALOGCXCM  I  )  )  (  1  .-XCA  (I)  )  *ALOGU  . -XCA  (  I  )  ) 
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ZB(I)=  XCB ( I ) *ALOG (XCBU))  f  Il.-XCBII))*  ALOG ( 1 .-XCB ( I ) ) 

HXS  (  I  )  =  XCR( I)  *  W 
GO  TO  270 

200  CON  =  -2.  +4.*  FMW 
205  YEMW  =  Y( I)  *  EMW 

A  *  CON  +  Y ( I )  -  YFMW 

B  =  2.  -  YEMW 

QUAD  *  B**2  +  4.*A«Y(I) 

IF  (QUAD)  107*210.210 

210  QUAD  =  SORT t QUAD) 

AA  »  2  .*A 

XCA{ I)  *  (-B  ♦  QUAD ) /AA 
IF  (XCA(I)-l. >220*230*230 
220  IF  IXCAU))  230.230.260 
230  ZB  (I )  =  XCAU  ) 

XCA(I)  =  (-R— QUAD) /AA 
IF  <  XC  A ( I ) )  250,255*240 
240  IF  (XCA(I)  -  1.)  260.250.250 
250  Z ( I )  =  XCA ( I ) 

255  Z A (  I )  *  2. 

GO  TO  280 

260  XCA1  *  1.  +  XCAU ) 

XCBU)  =  Y<I)/2.*XCA1/XCAU)  -  1. 

T01PX  =  2./XCA1 

TX01PX  =  T01PX  *  XCAU) 

ZA(I)=  TOlPX  *  (  XCAU  )*ALOG<  XCAU  ))  +  (  1. -XCAU  )  )*  ALOG(  1  .-XCA  (I )  )  J 
ZB(I)=TX01PX  *  {  XCBU)*AL0G(XC9U  )  )  +  (  l.-XCBl  I )  )*  ALOG(  l.-XCB  (I)  )  ) 
HXSU)  =  TX01PX  *  XCBU)  #  W 

270  Z (I )  =  ZA(I )  +  ZB (I) 

GXS< I )  =  HXS ( I )  +  Z  <  I > 

IF  (I  -  2499)280.  300  .300 

280  Y ( I  +  l )  =  Y{ I )  +  DY 

IF  IYII+11  -  YM*X)  290.290.298 
290  I  =  I  +  l 

GO  TO  { 105.205 ) .MODEL 

29F  YC I+l )  =  0. 

300  OO  380  J  =  1,1 

303  IF  (LINCT  -  42)  305.350.350 

350  LINCT  =  0 

WRITE  (6.360JMODFL  ,W 

360  FORMAT  (1H1.61X.  9HPROGRAM  Q  / /62X.  7HI42051)  / //44X .44HME2CY  PHA5 
IE  CALCULATIONS  -  EXCESS  QUANTITIES///  62X  ,  6HMODEL  .11/// 

2  57X.  3HW  =  E15.8  ////) 

305  IF  (XCBU))  330.310.330 
310  IF  ( Z A ( J ) -  1.)  320,315,320 

315  WRITE  (6,316)  Y(J)  ,ZB(J) 

316  FORMAT (Fll ,4,21X»31HNEGATIVE  VALUE  UNDER  RADICAL  =  E16.8.6X  , 

1  !5HI“AGINARY  POOTS  ) 

GO  TO 

320  IF  ( Z  A ( J ) -3 ,  )  325,327,225 

327  WRITE  <6.3?e)  Y(J)  ,  ZRU),  Z(J) 

328  FORMAT!  F11.4.13X,  11HXC3  ROOTS  (  ,E16.8  ,  3H  .  .E16.8  ,  40H  )  A 

IRE  NOT  IN  THE  RANGE  BETWEEN  0  AND  1  I 
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GO  TO  375 

323  IF  { ZA ( J )  -  2.)  330.323.330 

323  WRITE  (6.324)  Y(J)  *  ZBIJ)  .  Z(J> 

324  FORMAT  (F11.4.13X.  UHXCA  ROOTS  (  »E16.8  .  3H  .  .E16.6  *  40H  ) 

IRE  MOT  IN  THE  RANGE  BETWEEN  0  AND  1  ) 

GO  TO  373 

330  CALL  INT4D  (V.HXS.Y(J)  «TRP«  DHXS  ) 

CALL  INT4DIY.Z  ,Y( J) .TRP.DZ ) 

CALL  INT40IY.GXS.YIJ) .TRP.DGXS) 

370  WRITE  (6*371)  Y(J)  .XCA(J)  .  XCB(J)  .ZA( J J *ZB( J) *2 ( J)  .HXS(J). 

1  GXS(J).  DHXS  •  OZ  *  DGXS 

371  FORMAT  (1X.4HY  ■  E15.6.6H  XCA  *  E15.6.6H  XCB  »  E15.8.7H  ZA  • 

1  E15.8.7H  ZB  -  E15.8.6H  Z  *  E15.8/5X.6HHXS  ■  E15.B.8H  GXS  - 

2  E15.8.12H  DHXS/DY  -  E15.B.10H  DZ/DY  >  E15.8.12H  DGXS/DY  * 

3  E15.8) 

375  LINCT  ■  LINCT  ♦  2 
380  CONTINUF 


W  •  W+OW 

IF  (W>WMAXI  35.35.10 
ENO 

SENTRY  42051 


SDATA 


1 

•  96 

1.04 

.02 

.68315768 

.70315768 

.01 

2 

.96 

1.04 

.02 

.68315768 

.70315768 

.01 

2 

.05 

1.95 

.05 

.5 

9. 

.5 

1 

.05 

1.95 

.05 

.5 

9. 

.5 

1 

i 
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APPENDIX  ID. 


Computer  program  for  evaluating  the  high-temperature 
thermodynamic  properties.  For  theoretical  background,  refer 
to  AFML-TR-65-2,  Part  IV,  Volume  I  (1965) . (LenNole). 


JtaJJ 


$  IB JOB  GO.MAP 

SIBFTC  8064  LIST .REP.OECK 

C  EVALUATION  OF  HIGH  TEMPERATURE  THERMODYNAMIC  PROPERTIES 

C  VIA  NUMERICAL  INTEGRATION 

01  MENS  ION  DATA!  5001  tHEADI  12)  .TABLTdOO)  .KDIS«20)  »TABLY(  100)  * 

1  TABT ( 1 >  »TABY( 1 ) 

E0U1 VALENCE  I TABT .DATA! 1 ) )  ♦  ( TABY .DAT A( 201 ) ) .  t SST .DAT4I 401 >  I 
C 
C 

1  DO  2  I  >1.500 

2  DATA! I )■  0* 

CALL  SL1TEI0) 

CALL  AS138  (DATA. HEAD. NE) 

IF  (NE  -  216.5.4 

3  STOP 

4  WRITE  (6.5) 

5  FORMAT  (27H  (((INPUT  ERROR)))) 

GO  TO  1 

6  N  >  DATA (402 i 
DO  10  I  >1.200 

IF  ( TABT { I ) )4. 7 .10 

7  IF  ( TABT (1+1))4«11»10 

10  CONTINUE 

11  NOFT  «  1-1 
WRITE  (6.12)  HFAD 

12  FORMAT  (1H1.42X.27HAER0JEY  GENERAL  CORPORATION/  36X.41HLI QUID  ROCK 
1ET  PLANT. SACRAMENTO  CALIFORNIA  /  31 X.S5HEVALUAT ION  OF  HIGH  TEMPERA 
2TURE  THERMODYNAMIC  PROPERT 1 ES/ l 1H  DEPT.  4600.80X.8H JOB-8064/32X. 
312A6 ) 

WRITE  (6.13)  SST.  N.  (TABT ( I ) . TA3Y( I ) . TABT ( 1+50) »TABY( 1+50 ) » 
1TABT(I+100)»TABY(I+100)»TABTII+150>»TABY(I+1SOJ»  I  ■  1.50) 

13  FORMAT  (//.50X.20H  INPUT  DATA  //7X .3HSST .3X.7HSUB  INT/» 

1F10.2.8X.I2.//.23X.2HT  .8X.2HY  .10>*»10H  TI50I.10H  Y(50)»10 

2X.10H  T(100).10H  YdOO.lOX.lOH  T(150).10H  Y(150).//» 

3(15X.F10.2.F10.3.10X.E10.2.Fl0.3.10X,F10.2.ri0.3.10X.F10.2»F10.3 ) ) 

WRITE  (6.130) 

130  FORMAT  (1H1.4BX.21H  OUTPUT  DATA  //.10H  T  DEG  K  ,5X, 

15H  Y  .5X.5H  CP  »4X  .6HHT—H57 .4X.6HST-SST  »5X .5HF0GFE. 5X  » 1  OH  INTER 
2VALS ) 

KOlSd  )«  1 
I  ■  2 
K  ■  2 

TINTP  «  0. 

TABLTd)  «  TABT  ( 1 ) 

T  ABLY  1 1 )  ■  TABY(l) 

J  ■  2 

14  TABLTd)  •  TABT  ( J) 

TABLY(I)  -  TABY(J) 

IF  (TABTU+l)  )2P. 16.20 
16  IF  (TART( J+2) )4. 22.160 
160  J  *  J+2 
I  -  1  +  1 
KDIS(K)  «  J 
K  >  K  +  1 

TABLTd)  •  TABT(J) 

TABLY(I)  ■  TABY ( J) 

20  I  •  I  ♦  1 
J«J+4 
GO  TO  14 


64 


22  KDIS(K)  =  NOFT 
K  *  n 
LMN  *  1 

24  K  *  K  +  l 

NML  *  <0l$(K) 

LN  -  KOIS(K-H) 

IF  (LMN  -11231.231*25 
231  DELHT  •  0. 

DEL5T  *  0. 

CP  *  TABLY(l) 

FOGFE  »  DAT A( 401  ) 

GO  TO  27 

25  L«N  «  LMN+1 

DELHT  «  T ABLY ( LMN )  *  ( T ABLT ( LMN )  -  298*15) 

CALL  INT4D  ( TABT (NML > »TABY ( NML ) . TABLT C LMN) * YO.DY ) 

CP  »  DY*( TABLT (LMN)  -  298. 15 *  ♦  TABLY(LMN) 

TEMP  *  TABLY(LMN)*(1.  -  298.15/TABLT(LMN) ) 

CALL  SL ITET ( 1  * KL  ) 

IF  ( KL-1 >26*26.251 

251  CALL  INTGR  ( TABLT I LMN-1 ) .TABLT C LMN) *TX *N) 

CALL  INT4  ( TABT ( NML ).TABY( NML ).TX.YOFT) 

GOFTX  «  YOFT  *  (TX  -  298.15) /TX/TX 
CALL  INTGS  (GOFTX, T INT . .00001 *M) 

TINT  *  TINT  ♦  TINTP 
TINTP  «  TINT 

26  DEL5T  *  TEMP  +  TINT 

FOGFE  *  —DELHT /T ABLT ( LMN )+  DELST  +  S5T 

27  WRITE  (6.30)  T  ABLT  <  LMN ) »  T  ABLY ( LMN ) .  CP.  DELHT,  DELST.  fOGFE. 

30  FORMAT  (F10.2.F10.3.4F10.2.12X.I3) 

31  IF  (TABLT(LMN) )4,1 .32 

32  IF  ( T ABLT ( LMN )  -  TABT (LN ) ) 25. 34. 34 

34  IF  (TABT(LN)  -  TABT (NOFT ) ) 35. 1 .1 

35  CALL  SLITE( 1) 

GO  TO  24 

END 

5ENTRY  8064 

tOATA 
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APPENDIX  IV. 


Computer  program  for  calculating  the  free  energies  of 
the  monocarbide  phases  us  lag  the  Wagner  -Scholtky  vacancy 
model.  For  theoretical  background,  refer  to  AFM'^-TR-65-2, 
Part  IV,  Volume  1  (1965).  ..................  .  .(Jerry  Howard) 


*  *  O  u 


SEXECUTE 
IB JOB  807? 
IBFTC  807? 


IRJOB 

MAP 

LI^T.OFOC 


MONOCARBIDES  BY  USING  THE  VACANCY  MODEL, 

C8072-  PROGRAM  H.  THERMODYNAMIC  CALCULATIONS  FOR  GROUP  IV  METAL 
DIMFNSION  DATA (200) ,HtAD( 12 > *TEMP1 ( 20 ) »GA0(20 ) »*BOI 20 ) * 

1  ALGPAO ( 20 ) • ALGPBO ( 20 ) • GA ( 20 >  »GB( 20  > • RTLN ( 20 »  * 

2  ALFA ( 20) »TDGF ( 20 ) 

DIMENSION  A (20) *8(20) *BET (20) »GAM<  20) 

EQUIVALENCE  «CATA( 1) »TA) .(DATA (2 )»TAB> »(DATA(3) »TEL) » (DATA (4) ♦ 


1 

2 

3 

4 
8 
6 
7 


BETA) ,(DATA(5) »GAMMA), (DATA(6) »BETAP) , (DATA(7> .GAMMAP) , 
(DATA(E) *BETAPP) ,(DATA(9) .GAMMPP) ,(DATA) 11 
) ,XMt N ) , ( DATA ( 12 ) »XMAX ) » ( DATA (11) ,DELX  > . 

( DAT  A ( 1 5 ) * TEMPI )»{DATA(S5)»GA0) , 

( DAT A ( 7S ) .GEO ) , { DATA ( «*  >  , ALGPAO ) , ( DAT A ( 1 1 5 ) * ALGPBO ) , 

(  DAT  A  ( 1 5  5  )  ,  TICOi.  ;  )  ♦  ( DAT  A  (  38  )  t  AA )  ,  (  DAT  A  (  36  )  »  BB  )  *  (  DAT  A  (  37  ) 
»  A AP ) * ( DAT  A  ( 38 ) »9CP ) . (DATA (89) .AAPP ) • (DATA (40 ) ,3BPP ) 
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CALL  AS138 ( D AT A, HEAD t ME ) 
WR ITE( 6  *88 ) 

FORMAT (1H1) 

GO  TC  (3,4*2) »NE 
»'R  1  TC(  6 ,3  ) 

FORMAT ( ?0X,?4H  INPUT 
GO  TO  1 
CALL  EXIT 
PTAs.003 
ALN?  =  .6^1S 
R=l. 98726 
1*1 

ICNT1=T1C0NT+  ,01 


C  7  R  O  R  ) 


T  =  TEMPI ( I ) 
IFd.LT.TAE)  GO 
IF(T.LT.TBL)  GO 
BET ( I ) =9BTAPP 
GAM(  I  )»GA'*UPP 
Ad)  =  A^PP 
0(1)  s  Popp 

GO  TO  9 
BET ( I ) =BET  3 
GAM( I ) ■ GAMMA 
All)  =  AA 
B ( I)  =  BB 
GO  TO  9 


TC 

TC 
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— t- 


B€Td)-BET*P 
6 AN 1 1 l-GAMMAP 
Ad)  •  AAP 
Bill  «  BBP 
RT«R*T 

TOGF U  » -BET d ) *6 AM d ) *T 
GB  d  )  -A  d  )  *fi  d  )  *T 
R TLN ( I)  —GB d  ) -RT »ALN 2 
ALFA! I )-.5*EXPI-GB( 1 l/Rf) 

CAM!  o-RTt  N  d  )  -RT  *ALN2-T  DGF  ( I ) 


1-1*1 

IFd.LE.lCNTl) 


GO  TO  6 


20  J»1 

21  T-TEMP1IJ) 

CGA  »  GA(J) 

CGB  *  GBCJ) 

CRLNA  •  RTLN(J) 

CALF  -  ALFAIJ) 

RT  -  R*T 

WRITE (6 »22 )  HEAD 

22  FORMAT (1H1.42X.27HAEROJET  GENERAL  CORPORAT 1 0N/36X *4 1HLI QUl C  ROCKET 
1  PLANT .SACRAMENTO  CALIFORNIA//38X.12A6//19H  INPUT  DAT  A// > 

WRITE(6»100)  RET (J»  »GAM( J) *  A<J).  B(J) 

100  FORMAT (9X.16HL I NEAR  EQUATIONS/ 12X » 1 8H2  DELTA  GF  A  *  .E12.5. 

16H.  9  -  »E1 2»<?/i«?X»lflHGB*  A  *  .E12.5.&H.  B  «  .E12.5) 

WRITE ( 6 *101) GAO ( J) *GB0( J) .ALGPAOC J) .ALGPBO(J) 

101  FORMAT (9X»6HGA0  ■  .F7.0.8H.  GBO  ■  *F7.0»12H.  LOG  PAO  ■  »F7.3. 

U2H.  LOG  PRO  *  .F7.3) 

28  WRITE<6»29)  T 

29  FORMAT (9X.15HTEMPERATURE  IS  »F6*1»15H  DEGREES  KELVIN*/) 

WRITE! 6*102 )  CGA .CGB .CRLNA .CALF 

102  FORMAT <21H  OUTPUT  DAT  A//9X»$HGA+»  »E12.5»6H*6B*-  .E12.5* 
113H.RT  LN  ALFA-  *E12.5.7H.ALFA«  .E12.5//BX.1HX.11X.4HDGAB.10X. 
24HDGBB.11X.2HDG.12X.3HGAB.11X.3HGBB.12X.1HG.10X.6HLOG  PA.8X. 
36HLOG  PB) 

C1-4.576»T 
C2»C1*ALGPA0( J) 

C3*C1*ALGPR0<J) 

X-XMIN 


30  ARG  -  ABSIX/. 5-1.0) 
IF(ARG.LE«l.E-4)  GO  TO 
ir<X.LT.0.5)  GO  TO  34 


32 
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i 

i 


IFtX.GT.0.5)  GO  TO  35 

32  X  *  *300 

OCX  «  -CGA-RT»ALN?-CRLNA 
0GB  *  -CGB+CFA+DGA 
GO  TO  36 

34  DGA  ■  TDGF«  J)+CG3+RT*ALOG(  U.-?.*XJ/U.-X)) 

OGB  «  -CGB+RT*ALOG(X/n.-2**X)  I 

GO  TO  36 

35  06A  «  -CGA+RT*ALOG( (1.-X)/I2«*X-1.» > 

OGB  ■  TDGF(JI+CGA+RT*ALOG< <2.*X-1.»/X) 

36  VR1  -  1.  -X 
DG*VR1»DGA+0GB«X 
GAB*OGA+GAO( Jl 
GBB-DGR«-GBO(JJ 
G«VR1*GAo+X*G9R 
ALGPA-(DGA+C?)/C1 
ALGPB*(0GP+C31/C1 

WRITE(6.33I  X.0GA,DGR.0G.GAR.G«=R.G»AL6PA*ALGPR 

33  FORMAT (9C?X»F1?«5) ) 

X-X+OELX 

IF(X.LC.XMAX)  GO  TO  30 
J-J+l 

IF(J.LE.ICNTl)  GO  TO  21 

GO  TO  l 

END 
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